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£3 ■ Abstract 

We show that matrix Q x Q Self-dual type S-integrable Partial Differential Equations 
(PDEs) possess a family of lower-dimensional reductions represented by the matrix Q x 

ff") ■ uqQ quasilinear first order PDEs solved in [29] by the method of characteristics. In turn, 

these PDEs admit two types of available particular solutions: (a) explicit solutions and 

(y-j ' (b) solutions described implicitly by a system of non-differential equations. The later 

■ solutions, in particular, exhibit the wave profile breaking. Only first type of solutions 

is available for (l+l)-dimensional nonlinear S'-integrable PDEs. (l+l)-dimensional N- 

00 ■ wave equation, (2+1)- and (3+l)-dimensional Pohlmeyer equations are represented as 

examples. We also represent a new version of the dressing method which supplies both 
classical solutions and solutions with wave profile breaking to the above S'-integrable 
PDEs. 

^: 

.53 ; 1 Introduction 

There are several types of nonlinear equations of Mathematical Physics, which are referred 
to as integrable Partial Differential Equations (PDEs). We underline three following types of 
equations: 

1. Equations integrable by the Inverse Spectral Transform Method (ISTM), which are also 
called soliton or S'-integrable equations. Korteweg-de-Vries equation (KdV) has been 
discovered as the first representative of this class [1]. Among other methods of study of 
soliton equations we underline so-called dressing method, which was invented in [2, 3] and 
developed in set of papers, for instance, [4, 5, 6, 7, 8], see also [9, 10, 11]. 

2. Nonlinear PDEs linearizable by some direct substitution, or C-integrable equations [12, 
13, 14, 15, 16, 17]. Mostly remarkable is Hoph substitution allowing to integrate Burgers 
type nonlinear PDEs. 
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3. Quasilinear first order PDEs integrable by the method of characteristics [18] and by 
generalized hodograph method [19, 20, 21]. The later method is well applicable to (1+1)- 
dimensional nonlinear PDEs. Investigation of the higher dimensional systems by this 
method requires special efforts (see, for instance, the method of hydrodynamic reductions, 
which will be cited below). 

Many other method for study of the nonlinear PDEs have been developed. For instance, 
the algebraic-geometrical approach is well suited for the construction of local solutions [22, 23]. 
However, only restricted class of global solutions may be described in this way. A method 
for implicit description of solutions manifold is the method of hydrodynamic reductions [24, 
25, 26, 27, 28], where solutions of multidimensional integrable PDEs are described in terms of 
Riemann invariants of some (l+l)-dimensional first order quasilinear PDEs. The criterion for 
applicability of such technique to the given nonlinear PDE has been invented. This method 
allows to implement infinitely many arbitrary functions of single argument into solution of the 
nonlinear PDE satisfying this criterion. 

Our approach to the problem of construction of the particular solutions to the Self-dual type 
S-integrable PDEs is most similar to the method of hydrodynamic reductions. However, if it 
is applicable, we are able to describe the bigger solutions space. Namely, number of arguments 
in the arbitrary functions implemented into solution has no restriction and is defined by the 
particular nonlinear PDE. But we have no criterion for applicability of this algorithm to the 
given nonlinear PDE, unlike the method of hydrodynamic reductions. 

We exhibit a class of solutions with wave profile breaking to the S'-integrable PDEs selected 
above. This phenomenon is not surprising, since, simultaneously, these solutions are solutions 
of the well established class of first order quasilinear PDEs in lower dimension, which are known 
to possess solutions spaces implicitly described by systems of non-differential equations [29, 30]. 
Remark that solutions with breaking of wave profile to dispersion-less Kadomtsev-Pertviashvili 
equation (dKP) (which is compatibility condition of appropriate pair of vector fields) have been 
obtained and described in [31] using ISTM. 

In the rest of Introduction we, first of all, recall the method of characteristics applicable to 
a certain type of matrix first order quasilinear PDEs, Sec. 1.1. Then we remind some important 
properties of the Self-dual type nonlinear S'-integrable PDEs, Sec. 1.2. Basic novalties appearing 
in this paper are underlined in Sec. 1.3. Finally, we represent the brief contents of the paper in 



1.1 First order quasilinear PDEs integrable by the method of char- 
acteristics 

As it was mentioned above, we associate the given Self-dual type S-integrable nonlinear PDE 
with the family of lower dimensional first order matrix PDEs integrable by the method of 
characteristics. Let us describe this family of PDEs. Throughout this paper we will write 
superscripts in parentheses in order to distinguish them from powers. 
It is well known that the scalar first order PDE 



can be solved by the method of characteristics [18]. For instance, solution u(x, t) of (1) for 



Sec.1.4. 



N 




(1) 



fc=l 
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p = is defined implicitly by the non-differential equation 

u = f( Xl -pW(u)t,..,x N -pW(u)t). (2) 

Vector generalizations of equation (1), i.e. the systems of several coupled first order scalar 
PDEs, have been investigated in set of papers [19, 21, 32, 33] using generalized hodograph 
method in most cases. Recently [29] the matrix generalization of eq.(l) has been solved using 
algebraic approach: 

N 

wt + J2 w * k P (k) ( w ) = P( w ) + K T PM], (3) 
fc=l 

where w is the unknown QxQ matrix function of the N+l independent variables (xi, t) G 
R^" 1 " 1 , T is any constant diagonal matrix, [•, •] is the usual commutator between matrices, and 
p^ k \ p, p : R ^ R, k — 1, .., N are N + 2 arbitrary scalar functions representable by the 
positive power series: 

i>0 i>0 i>0 

where and are scalar constants, so that the quantities p^ k \w), k = 1, .., N, p(w) 

and p(w) are well-defined functions of the matrix w. 

Remark. We may always apply the following change of variable t and field w 

N 

d t + J2 -ft, ™ - e-^ )Tt we &i0)Tt (5) 

3=1 

which is equivalent to putting 

a (j;0) = Q) ~(o) = (6) 

in the eqs.(4) without loss of generality. 

Solutions space to the matrix equation (3) with p = is implicitly described by the following 
algebraic equation: 

W a/ 3 = (7) 
Q 

(e-r^F^xJ - p {1 \w)t, ...,x N I- pW{w)t)e T ^ 



7=1 

a — 1, ... ,Q, j3 — 1, . . . , Q, I is Q x Q identity matrix. 
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If p = 0, then eq.(3) has infinitely many commuting flows. In this case, it is convenient to 
replace eq.(3) by the next equation 



N 



^ + E^^ ) H = K^ H H1. m=l,2,..., (8) 



fc=i 



where index m enumerates commuting flows, p( mfc ) are arbitrary scalar functions representable 
by a positive power series 

p^\z) =J2a (mr ' l) z\ p {m \z) =J2® {m ' i)zi , (9) 

i>0 i>0 
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a ( m r^) ^ q,("»;») are scalar constants, are diagonal matrices. Note, that the transposed 

equation 

N 

w tm + J2p {mk \™)ti Xk = \p {m \w)T m ,w], w = w T , m = l,2,... (10) 
k=i 

is also integrable. 

Equation (3) as well as its generalizations whose solutions spaces may be completely de- 
scribed implicitly by the system of non-differential equations will be referred to as PDEO in 
this paper for the sake of brevity. 

Using either the algebraic manipulations (Sec. 2) or the modified version of the dressing 
method (Sec. 3), we will show that PDEOs (8) may be viewed as lower dimensional reductions 
of the classical Self-dual type S'-integrable PDEs. These PDEs are referred to as PDEls in this 
paper. 

1.2 Family of PDEls 

PDEls may be viewed as compatibility conditions of the following overdetermined system of 
linear PDEs for some spectral function V(X;x) (spectral system): 

VtjX] x) + Yl Xnm ' V ^ X )+Y1 ^ m *V(\; x)T® = AV(A; x)qW(x), (11) 

i=l i=l i=0 

m = 1,2. 

Here all functions are Q x Q matrix functions, A is a complex spectral parameter, n mi and 
h mi are some integers, N^X = max(n mil ,n mi2 , i x = 1, . . . , N, i 2 = 1, . . . , N); g (mi) (x) are 
some matrix functions of x — (xi, x 2 ..., ti, t 2 ...). Throughout this paper we use Greek letters in 
the lists of arguments for spectral parameters. The feature of this system is that its LHS is a 
first order differential operator having coefficients which are polynomial in A and independent 
on the vector parameter x, i.e. x-dependent potentials appear only in the RHS, which is 
also polynomial in A. Compatibility condition of any pair of eqs.(ll) results in PDE1 for 
the functions q( mi \ Note that eq.(ll) uses only multiplication by the spectral parameter A and 
does not involve derivatives of the spectral function with respect to the spectral parameter. For 
this reason such S-integrable PDEs as dispersion-less Kadomtsev-Petviashvili equation (dKP) 
and Heavenly equation may not be considered in the framework of the represented dressing 
algorithm. 

Remark, that the spectral system (11) may be replaced by the equivalent one 

N N 

d tm V{\- x) + Y, A n ™ (A, u) * d Xi V(u; x) + ^ A^(X, u) * V(v; = (12) 

i=i i=i 

N (rn) 
1 v max 

A l {X,u)*V{u;x)q (mi \x), 

i=0 

where function A(X, /i) replaces multiplication by A and we define operator A % as follows: 

A 1 = A* ■ ■ ■ * A. (13) 
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Emphasise that A is independent on x. The eq.(12) reduces to the eq.(ll) if A(X, v) = X5(X — 
fj,)I, where I is Q x Q identity matrix. 



Example 1: iV-wave system. Linear spectral system reads in the simplest case: 

d tm V(X;x)- A(X,u)*V(iy;x)T {m) = V(X;x)[T^ m \w (0 \x)}, m = l,2. (14) 
Compatibility condition of this system yields the next PDE1: 

prw^g)] _ [T w,wV>] + [[r< 2 >,«,< 0) ], [rw «,<">]] = o. (is) 

Most physical applications require reduction w^p = wp°a, where a ^ (3 and "bar" means 
complex conjugate. This equation is well known as (l+l)-dimensional iV-wave equation [10]. 
We consider eq.(15) without reductions. 

Example 2: Pohlmeyer equation. Another example is assotiated with the following spec- 
tral system: 

d tm V(X; x) + A(X, v) * d Xm V(u; x) = V(X; x)w^(x), m = 1,2. (16) 
Compatibility condition of this system yields the next PDE1: 

i-A^M, ( 17 ) 

which may be written in a different form: 

(J-'Jt^-iJ-'Jtzh^O, (18) 

This is Pohlmeyer equation [34, 35, 36]. Most applications in Physics is assotiated with its 
reduction J + = ±J, det J = 1, which yields (Anty) Self dual Yang-Mills equation (ASDYM). 
This equation has been studied by many authors: [37, 38, 39, 40, 41, 42, 43, 44] (instanton 
solutions), [45] (merons), [46, 47] (finite gap solutions). (2+l)-dimensional version of ASDYM 
assotiated with the reduction 

d X2 J = d h J, d X2 w^ = d tl w^ (19) 

is also well known: 

«fi - wgi, = [w£\wj®] or (20) 
(J Jti)ti ~ {J Jt 2 )x 1 — 0. 

It has been studied in [48, 49] (Initial Value Problem), [50, 51] (localized solutions). Here we 
will deal with eqs.(17) and (20a) without any reduction for the field w^°\ 
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1.3 Basic novalties 



There are two manifolds of particular solutions to PDEls. First manifold is assotiated with the 
uniquely solvable linear integral equation for some spectral function (Zakharov-Shabat method 
[52], classical 9-problem [53]). We concentrate on the second manifold of solutions which is 
described implicitly by a system of non-differential equations (see Sec. 2). Simultaneously, this 
manifold is solution manifold to appropriate lower dimensional PDEO. Proper choice of initial 
data leads to wave profile breaking. Such solutions for Pohlmeyer equation will be discussed. 
In turn, second manifold of solutions has a sub-manifold of explicit solutions. Regarding (1 + 1)- 
dimensional iV-wave equation, only such solutions from the second manifold are available. It is 
likely that the same conclusion is valid for any other (l+l)-dimensional S'-integrable model. 

It is remarkable, that there is a version of the dressing method, which joins both manifolds of 
solutions to DPEls, Sec. 3. This fact is remarkable, because the second manifold is beyond the 
scope of the classical dressing method. In addition, this fact demonstrates that two significantly 
different classes of PDEs (S'-integrable PDEls and integrable by the method of characteristics 
PDEOs) are combined by the dressing method. There are examples of similar combinations. For 
instance, dressing method combining S- and C-integrable models has been suggested in [54]; 
dressing method combining C-integrable equations and equations integrable by the method of 
characteristics has been proposed in [30]. 

All these examples tell us about flexibility of the dressing method which is important for 
development of the uniform method for integration of multidimensional nonlinear PDEs. 

Usually, we will use notations PDEl(t mi , t m2 ; it/ )), PDE0(t m ; no, w) and PDE0(t m ; ko, w) 
with arguments reflecting field as well as derivative(s) of field with respect to variable(s) t rn 
appearing in PDE. Of course, PDEOs and PDEls involve also derivatives with respect to x^. 
However, we do not represent them in the lists of arguments. Here w(x) is n Q x n Q matrix 
and w(x) is koQ x koQ matrix block functions: 



w(x) 



w<®(x) 
I 




w(x) = 



w 



(1), 



I 



_yj(k -2) { 



x) 





-w<®(x) 



X 





I 







w 



(n -2), 






-w 



x) 








w { - no ~ 1 \x) 
r (i) 

r (2) 
r (no-l) 



(21) 



where / and are Q x Q identity and zero matrices, and it)W are Q x Q matrix functions, 
n and k are arbitrary integers, i = l,...,n — 1 are arbitrary constant Q x Q 

matrices (if = in eqs. (25,26)) or arbitrary diagonal Q x Q matrices (if ^ 0). Each 
PDEl(ti, ti\ w^) written for the matrix field it/ -* and involving derivatives with respect to t\ 
and t2 possesses a family of PDE0(t m ; no, if) and PDE0(t m ; ko,w), m = 1, 2, no, ko — 1, 2, . . . , 
as lower dimensional reductions. Of course, PDE0(£i; n ,w) is compatible with PDE0(t 2 ; n ,w), 
as well as PDE0(£i; no,w) is compatible with PDE0(t 2 ; n ,w). 

For instance, in Sec. 2 we will derive the following PDEOs assotiated with iV-wave equation 
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(Sec.2.1.1): 

w tm -[w,TW W } = 0, (22) 
tit m -[wf£\w} = 0, m=l,2, 

T^ m) = diag( T (m) - - -T (m) j , f( m) = -Ti m) 

n 

and PDEOs assotiated with Pohlmeyer equation (Sec. 2. 1.2): 

w tm + w Xm w = 0, (23) 
w tm +ww Xm = 0, m = 1,2. 

If reduction (19) is imposed, then eq.(23) reduces to 

iu tl + w xi w = 0, w t2 + w tl w = 0, (24) 
W tl + = 0, w t2 + ww tl = 0. 

General forms of PDE0(t m ; no, w) and PDE0(t m ; ko, w), m — 1, 2, which are lower dimensional 
reductions of some PDEl(ti,t 2 ;iu (0) ) are following: 

% + E^./ )(m!) H = Ei w ' T ^ W N. m = l,2, (25) 

i=l i=l 
N N 

w tm + Y.P {m ^^)^-Y}P { ^^) f kl^l m = l,2, (26) 
i=i i=i 

where functions p(™) are defined by the series similar to the eq.(9): 



^\z) = ^a {mj ' l) z\ (27) 



P 

i>0 



and a {mj '^ are scalar constants. Remarks on the origin of eqs.(25) and (26) will be given in 
Sec. 3. 3, see n.3,5 and n.4,6 respectively. Eq.(26) has the form of transposed equation (25) with 
replacements 

w = w T , T^^fW. (28) 

Because of this similarity, we will deal with PDEl(ti,t2',w^) and PDE0(t m ; n , w), m = 1,2. 
Regarding PDE0(t m ; k ,w), only some details will be given. 

Remark 1. Although any PDE1 is compatibility condition of eqs.(12), we are not able to 
represent general form of PDEls explicitly, unlike general form of PDE0, eqs. (25,26). 

Remark 2. The eq.(25) defers from the eq.(8) by the sum in the RHS. However, eq.(25) may 
be treated by the methods developed in [29, 30] for eq.(8). 

Remark 3. The eqs. (22,23,25,26) (where m — 1,2) have infinitely many commuting flows 
corresponding to m > 2 in these equations. 
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1.4 Brief contents 



In the next section (Sec. 2) we give algebraic description of PDEl(ti, t 2 ; w^ '), PDE0(t m ; no, w) 
and PDE0(t m ; k ,w). Derivation of PDE1 (namely, (l+l)-dimensional iV-wave and (3+1)- and 
(2+1)- dimensional Pohlmeyer equations) assotiated with appropriate PDE0(t m ; n , w) and 
PDE0(t m ; ho, w), m — 1,2 is given in Sec. 2.1. Assotiated linear overdetermined systems have 
been described therein as well. Using slightly modified results of [29, 30] we represent the non- 
differential matrix equations implicitly describing the solutions spaces to the above PDEOs and 
assotiated solutions manifolds to the appropriate PDEls, Sec. 2. 2. We derive some manifolds of 
explicit solutions to iV-wave and Pohlmeyer equations and describe solutions with wave profile 
breaking for (2+1)- and (3+l)-dimensional Pohlmeyer equations. 

A new version of the dressing method describing PDEOs and their relations with PDEls is 
proposed in Sec. 3. We have derived examples of PDEOs and examples of PDEls (following Sec. 2, 
iV-wave and Pohlmeyer equations are taken as examples of PDEls) together with appropriate 
linear overdetermined systems (spectral systems), Sec.3.1. Classical solutions manifolds to 
PDEls (Sec. 3. 2.1) and solutions manifolds with wave profile breaking to PDEls and PDEOs 
(Sec. 3. 2. 2) have been derived in Sec. 3. 2. 

Although this version of the dressing method is mostly straightforward in spirit of the 
derivation of PDEls, it does not supply full solutions spaces to PDEOs, unlike the method of 
characteristics. The second version of the dressing algorithm exhibits fullness of the solutions 
spaces to PDEOs, Sec. 3. 3. 

Conclusions are given in Sec. 4. 



2 Algebraic description of PDEls and PDEOs 

We consider three examples of PDEl(ti, ti\ w^) (namely, eqs.(15) and (17) together with its 
(2+l)-dimensional reduction (20a) ) and appropriate families of PDE0(t m ; n , w) and PDE0(t m ; k , w) 
m = 1, 2, in the framework of the algebraic approach, see Sees. 2. 1.1, 2. 1.2. The solutions spaces 
to PDEl(ti,i 2 ;iu (0) ) and PDE0(£ m ; n , w), m = 1,2, will be investigated in Sec.2.2. 



2.1 Derivation of PDEls and appropriate families of PDEOs 
2.1.1 (l+l)-dimensional iV-wave equation 

It is well known, that any S-integrable equation can be represented as system of two commuting 
discrete chains. Mostly explicitly this representation is given in Sato approach to integrability, 
see for instance, [55]. However, any known dressing method gives rise to this representation. 
We show, that iV-wave equation can be derived from the following pair of commuting discrete 
chains with two discrete parameters: 

W ^ = w (Hn+l)) T ( m ) _ T (m) w ((k+l)n) + w (fc0) T (m) w (0n) ? m = 1? 2 , (29) 

supplemented by the next non-differential relation among w^: 

w ((k+l)n) = w (k(n+l)) + w (k0) w (0n) ^0) 

(see also Sec. 3. 1.3, eqs. (150,151) with = 0). In fact, first of all remark, that eq.(29) gives 
rise to two alternative discrete chains with single discrete parameter in view of (30). The first 
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chain follows after putting k — and eliminating -u/ ln ) using eq.(30): 

wt] -[w {n+l \T^] + [T {m \w i - Q) ]w { - n) = m = l,2, (31) 

where 

w {n) =w (0n \ n = 0,l,.... (32) 

The second discrete chain follows after putting n = and eliminating from (29) using 

(30): 



where 



- [* (fe+1) ,^ (m) ] +* (fc) [* (0) > T(m) ] = 0, m=l,2, (33) 



(*) = «,(*»), A; = 0,1,..., ^°)=^°). (34) 



Finally, in order to write PDE1 for the function w^°\ we fix n = in (31) and eliminate 

or fix k = in (33) and eliminate In both cases we will end up with the same equation 

(15). 

Along with PDEl(*i, t 2 ; w (0) ) (15) we may derive the family of PDE0(t m ; no, w), m — 1, 2, 
from the eqs.(31) imposing the reduction 

n - l 

iu^x) = w {i) (x)r^, (35) 

i=0 

where no is an arbitrary integer and are arbitrary diagonal constant matrices. Similarly, we 
may derive PDE0(t m ; k ,w), m = 1,2, from the eq.(33) with the reduction 

fco-l 

w^ k °\x) = f Xt) w (t) (x), w (0) = w (0) , (36) 

where ko is arbitrary integer and are arbitrary diagonal constant matrices. Both (35) 
and (36) are closures of the chains (31) and (33) respectively. However, parameters and 
f(°) may be removed from the PDEOs by the shifts of fields ■u/™ 0-1 ) + — > u^™ 0-1 ) and 
—yj(ko-i) + f(o) _^(fco-i)_ Thus hereafter we take 

r<°> = f <°> = (37) 

without loss of generality. Parameters no, ko, and f ^ appear in the definitions of w and w, 
eqs.(21). The simplest explicit examples of PDE0(t m ; n , w) and PDE0(t m ; k , w) are following: 

PDE0(t m ,l,w): wl°2 + [T^ m \w^}w^ = 0, m = l,2, (38) 
PDE0(t m] 2, u>), r« = : wf} - [w (1 \T (m ^] + [T<- m \w<®] W W = 0, (39) 

™ [ t2 + [T (m) , w (0 V (1) = 0, m = 1, 2, 
PDM)(t m ; 1,«>) : «>£! + * (0) = 0, (40) 
PDE0(t m ;2,w), f« = 0: - T^] + w^ ), T^] = 0, 

4l ) +* (1) [w (0) ,T (m) ] =0, m = l,2. 



It is not difficult to observe that the above PDE0(t m ; n , w) and PDE0(t m ; k , w) with arbitrary- 
no and k , may be written in the form (22a) and (22b) respectively, where w and w are defined 
by the formulae (21). 

Since the eq. (22b) follows from the eq.(22a) after transposition with replacements (28), it 
is enough to study one of them, say, eq.(22a). Slightly modifying result of [30], we describe the 
solutions space to the eq. (22a) with arbitrary integer no by the next implicit algebraic equation: 

^ ( - t {T^) aW t m E(T^Vt„A 

w aP = 2_^(e m=i F aj (w) e-=i J , (41) 

7=1 

a,/3 = 1, . . . ,n Q, 

where F(z ) is n Q x n Q matrix function of single argument. Details of derivation of this 
formula by the dressing method are given in Sees. 3. 2. 2 and 3.3. However, the structure of w 
tells us that F must have the following structure: 

_ , s f arbitrary scalar function of arguments, a < Q , . nS 

^ = I 4**, a > Q • < 42 » 

so that equation (41) becomes an identity for a > Q. Thus, the square matrix algebraic equation 
(41) reduces to the rectangular matrix equation with a = 1, . . . , Q and (3 = 1, . . . , rioQ. The 
scalar functions of single argument F a/3 (z ), a = 1, . . . , Q, (5 = 1,..., n Q are arbitrary. The 
simplest case n = 1 yields w = -u/°\ = T^: 



w «l = E ( e F ^(^ (0) ) )^ a, p = 1, . . . , Q. (43) 



7=1 



Assotiated linear overdetermined system of PDEs. As it was mentioned above, the 
linear spectral problem for iV-wave equation is given by the eqs.(14), where V(X; x) is, generally 
speaking, a rectangular IQ x Q matrix function, I is some integer (to anticipate, / is assotiated 
with the dimension of the kernel of the integral operator in eq.(114); / = 2 in Sees. 3. 1,3. 2). One 
can extend this spectral system introducing discrete parameter n by the following formulae: 

A{X,u)*V {n) {u;x) = V (n+1) {X;x) + V {0 \X;x)w {n \x), (44) 
V t ^\X;x)-A(X,u)*V {n \u;x)T^ m) = \/ (0) (A; x)[T {rn \ w {n \x)}, (45) 
n = 0,l,2,...,m = 1,2. 

which becomes eq.(14) if n — and V = V^°\ This extension can be derived formally, for 

instance, using a version of the dressing method developed in Sec. 3, see Sec. 3. 1.3 eq.(149) with 

s (m) = Q 

Let us consider the reduction (35,37) for the fields w^ n \ which causes the appropriate re- 
duction for the spectral functions : 

"o- l 

y(n ) = ^ V {i) r {i) . (46) 
i=i 

In view of this reduction, we may introduce the block-vector spectral function 

V= [y (0) .--V (no ~ 1) ]. (47) 
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Then eqs. (44,45) can be written in a compact form: 



A(\,v)*V(v;x)=V(\;x)w(x) (48) 
V tm (A; x) - A(X, v) * V(i/; x)T^ = V(A; x) [T^\ w(x)}, m = 1, 2. (49) 

An important feature of this system is the eq.(48) which has no derivatives with respect to 
variables ti and X{. Thus, w[x) in the RHS of this equation may be treated as a matrix 
parameter. For this reason, eq.(48) may be explicitly solved, at least, for diagonalizable matrix 
w (more complex case, when w is transformable to the Jordan form, is not considered here): 

w(x) = P(x)E(x)p- 1 (x), (50) 

Paa = 1, (51) 

where E is a diagonal matrix of eigenvalues, P is a matrix of eigenvectors normalized by the 
condition (51). 

Let A be in the following form: 

A(X, v) — X8(X - v)I h (52) 
where // is IQ x IQ identity matrix. Then solution of eq.(48) reads: 

V(X;x)P(x) =V(x)5(XI - E), (53) 

where 

V = [V {0) ■ ■ ■ V( n °-V] (54) 

Thus, V(A; x)P(x) is delta-function of A. This is an interesting fact and a non-trivial reduction 
imposed on the spectral function V(A;x). It is remarkable, that there is a dressing algorithm 
providing consistency of this reduction with eq.(49), see Sec. 3. 1.3. 

There is an important remark. Remark: The formula (53) means that 

V(A;x)^0 as A^oo (55) 

for any function w(x). Such a behaviour of the spectral function is in contradiction with the 
classical dressing, where the spectral function is not vanishing at infinity. For this reason the 
solutions described in this section are not observable in the framework of ISTM. 



2.1.2 Pohlmeyer equation 

One can show that Pohlmeyer equation may be derived from the following pair of commuting 
discrete chains with two discrete parameters: 

wf n) + w« k +1)n) - w^wW = 0, (56) 

supplemented by the eq.(30) (see also Sec.3.1.3, eq.(150,151) with = 1, = 0). First 

of all, similar to the Sec. 2. 1.1, eq.(56) gives rise to two alternative chains with single discreet 
parameter. The first chain appears after putting k — and eliminating -u/ 1 ™) using eq.(30): 

™& ) +™£r 1) +™& (n) = > ™ = 0,1,2,..., m = l,2, (57) 



11 



where fields are given by the eq.(32). The second discrete chain appears after putting 
n = and eliminating using eq.(30): 

+ - * W ^°2 = 0. n = 0,1,2,..., m = l,2, (58) 

where fields are defined by the eq.(34). Finally, in order to write PDEl for w^°\ we put 
n = in eq.(57) and eliminate or put k = in eq.(58) and eliminate In both cases 

we end up with eq.(17). 

Along with PDEl(ti, t 2 ; u> (0) ) (17) we may derive PDE0(t m ; no, w), m = 1,2, imposing the 
reduction (35,37) to the eq.(57) or PDE0(t m ; ho, w), m— 1,2, imposing the reduction (36,37) 
to the eq.(58). The simplest examples are following: 

PDE0(t m ;l,w): w^+w^w^=0, m = 1,2, (59) 

PDE0(t m ; 2, w), r« = : + + w^w™ = 0, 

+w^w (1) =0, m = 1,2, 

PDE0(t m ,l,w): w<®-w®wW=0, m = l,2, (60) 

PDEQ(t m ; 2, w), = : wf} + w<& - w^w^ = 0, 

wf ] - w {l) wf ] =0, m = 1,2. 

It is not difficult to observe that the above PDE0(t m ; n , w) and PDE0(t m ; /c , w) with arbitrary 
n and k may be written in the form (23a) and (23b) respectively, where w and w are defined 
by the formulae (21). 

Similar to the previous section, we represent the detailed consideration only to the eq.(23a). 
Its solutions space is described by the next non-differential equation [29]: 

= y2[F ai (w,x 1 I no -wt 1 ,x 2 I no -wt 2 )) , (61) 

7=1 

a, (3 = 1, . . . ,n Q. 

Details of derivation of this formula by the dressing method are given in Sees. 3. 2. 2 and 3.3. 
Here, similar to the iV-wave equation, F must have the structure predicted by the structure of 
w given by the eq.(21a): 

_ , x f arbitrary scalar function of arguments, a < Q 

Fa ^ Z ^ Z2) = \6 a ,z , a>Q ' (62) 

i.e. the square matrix equation (61) reduces to the rectangular one with a = 1, . . .,Q, (3 = 
1, . . . , n Q and arbitrary scalar functions of three arguments F a/3 (z , z 1 ,z 2 ). The simplest case 
no — 1 yields w = 

Q 

S = E (^(^ ^i' " ™ (0) ti, x 2 I - w (0 H 2 ^ , a, (3 = 1, . . . , Q. (63) 



w 

IXfJ 

7=1 



Let reduction (19) be used. Consider PDE0(t m ; n , w) (24), m = 1,2, corresponding to 
PDEl(ti,i 2 ;«;W) (20). Its solutions space is implicitly described by the next algebraic equation 



n Q 



a/3 = 52 [F ai (w,xiI no - Wti + w 2 t 2 ) ) , (64) 

, V / -y/3 

7=1 

a = 1, ... ,Q, (3 = 1,.. .,n Q 
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with arbitrary scalar functions of two arguments F a p(zo,z-\). The simplest case no = 1 yields 
(w = w<®): 

Q 

w S = E( F ^( w(0) ' Xl/ - w(0)tl + ( w(0) ) 2t2 )) ' «,/3 = l,...,Q. (65) 

7=1 7/3 

Assotiated linear overdetermined system. The classical overdetermined system for Pohlmeyer 
equation is given by the eq.(16), where V is IQ x Q matrix function, / is some integer. System 
(16) may be extended as follows (such extension may be derived, for instance, using a version 
of the dressing method developed in Sec. 3, see Sec.3.1.3 eq.(149) with = 1, = 0): 

V£\\;x) +A(\,u) *V<g(v;x) - V<°> (A; *)«;£>(*) = 0, (66) 

771= 1,2, 

where \/ (n+1) is related with \/ (n) by the eq.(44). Eq.(66) becomes eq.(16) if n = and V ee V^. 

We consider the reduction (35,37,46) and introduce the block- vector spectral function (47), 
which allows us to write the eq.(66) in a compact form: 

V im (A; x) + A{\, v) * V Xm (u; x) - V(A; x)w x Jx) =0, m = 1, 2, (67) 

while eq.(44) becomes eq.(48). Thus, the spectral equation (48) appears in both spectral 
systems corresponding to iV-wave and Pohlmeyer equations. This spectral equation is universal 
equation assotiated with reduction (46). Suppose that the matrix w is diagonalizable and 
representable in the form (50). Let A(X, v) be in the form (52). Then V(A; x) is given by the 
eq.(53). The dressing algorithm consistent with reductions (35,37,46) will be represented in 
Sees. 3. 1.3 and 3.3. 



2.2 Solutions space 

As before, we assume diagonalizability of w, i.e. w may be represented in the form (50), and 
describe two different types of solutions to PDE0(t m ; n , w), m — 1, 2, and assotiated solutions 
to PDEl(*i,f 2 ;iu (0) ): 

1. All eigenvalues of w are independent on x. This leads to the explicit solutions w without 
the wave profile breaking. In particular, the solutions in the form of rational function of 
exponential functions with linear in x arguments are available. Solitary waves belong to 
this type. 

2. Some or all eigenvalues of w depend on x. This leads to the solutions w with wave profile 
breaking. Such specific behaviour is typical for the systems of hydrodynamic type. 

In this section we will use the description of the solutions space to PDEOs given in [29], i.e., 
instead of the equations (41,61,64) we will consider the appropriate non-differential equa- 
tions describing the solutions spaces to the first order quasilinear PDEs for E and P. For 
PDE0(t m ; no, w) written in general form (25) these quasilinear PDEs read: 

N 

E tm + J2E Xi p {mi \E) = 0, m = l,2, (68) 

i=l 

N N 

Pt m +Y, P **P imi) ( E ) = Y,^ T n^P {m HE), m = l,2. (69) 
i=i i=i 
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This system is solvable by the method of characteristics [29] giving the following non-differential 
equations for Ep and P a p: 

Ep = (E^x 1 -^p^ 1 \E^t m ,...,x N -^p^ N \E p )t m Y (70) 

\ m=l m=l / 

2 N 

Pa/3 = e ™=^ =1 i Ep, xi — ^ p^ ml \Ep)t m , (71) 



m=l 



2 AT 



^-^P (mJV) W m ) e m =^ =1 , a,(3 = l,...,n Q, 



m=l 



where F^ and F^ are arbitrary functions of arguments if only there is no restriction on the 
structure of w. However, unlike [29, 30], we have the particular (rather then general) matrix 
structure of w, which requires the appropriate particular matrix structure of P. In fact, eq.(50) 
may be written in the form 



n l. 



^2 w ai P lg = PapEp, a, = 1, . . . , n Q. (72) 

7=1 

Eqs.(72) with a = 1, . . .,Q and (3 = 1, . . . ,noQ should be taken as definitions of j = 
0, . . . ,n — 1 in terms of E and P. Other equations of the system (72) yield the next linear 
relations among the elements of P: 

noQ noQ 

Y,w ai P lP = P a pEp P {a - Q) p + r< &y 1)p iP = p *P E Pi (73) 

7=1 7=(n -l)Q 



, a — a — sQ, a — Q + 1, . . . ,n Q, (3 = 1, . . . , n Q, 



where [•] means an integer part of number. Eq.(73) with Ep ^ and arbitrary may be 
always solved for the following set of elements P a p\ 

V = {P(a-Q)?, Pa(a-Q)] « = Q + 1, . . . , n Q , (3 = 1, . . . , n Q, (3 + a - Q}. (74) 

Other n\Q 2 — (no — l)n Q 2 — n Q = rioQ(Q — 1) elements of P are defined by the eq.(71) 
with arbitrary scalar functions F^g ■ The case uq > 1, Ep = for some (3q results in some 
restrictions for so that solution to the eqs.(73) will differ from V. This case will not be 
considered. 

All in all, the following algorithm gives either the explicit solutions to PDE0(t m ; n , w) and 
assotiated PDEl(£i, t 2 ; w^) (E(x) = const) or the non-differential equations (E(x) ^ const) 
describing the solutions spaces to PDE0(t m ; n , w), m — 1,2, and assotiated PDEl(ii, t 2 ; w^): 

1. Fix the arbitrary eigenvalues Ep ^ 0, a — 1, . . . , n Q, in the form (70). 

2. Fix n Q(Q — 1) elements of P a p V in the form (71). 



14 



3. Find other elements of P G V solving the linear system (73) with normalization (51). 

4. Find w using formula (50). 

We will describe some solutions to iV-wave and Pohlmeyer equations in Sees. (2.2.1) and (2.2.2) 
respectively. 

Implicit representation of the solutions space by eqs. (70,71) reveals solutions with wave 
profile breaking. This phenomenon is assotiated with Ep(x) ^ const (Pohlmeyer equation). 
We refer to the break points manifold related with Ep as Ep. We will show that break points 
manifold for P a p is also Ep. Then the break points manifold corresponding to w will be £ = 



I \noQc 



(75) 



2.2.1 iV-wave equation 

The equations (68,69) corresponding to the eq.(22a) read: 

E tm = 0, 
P tm = [P,T^}E, m=l,2. 

Then eq.(41) describing the solutions space to iV-wave equation will be replaced by the next 
pair of equations (the appropriate variant of the eqs. (70,71)): 

E a = e a = const, a — 1, ... , noQ, 



(76) 



P a p = e -=i (io)a/3e m=1 



a 



1,...,Q, p = l,...,n Q 

where P 0a p are constants, related by the eq.(73). Only n Q(Q — 1) of them may be arbitrary. 
Appropriate function will be rational function of exponents. 

An important problem is to exhibit such solutions which have no singularities in domains 
of all independent variables. Two examples are given below: 

1. Let e a = ie a , e a be real parameters, and require detP ^ 0, t± G R, t 2 G K. For instance 
let Q — 3, n — 1, P 23 — -P032 — 0, other elements of P equal i, i 2 = —1 (Po a /3 — i), 



= a n , e a = i (a — 1). Corresponding P reads: 



1 







1 



(77) 



which produces periodic w 



(0) 



2. Let e a be real parameters, detP 7^ 0, t 1 G R, t 2 G R. However, unlike the previous 
example, this is not enough to obtain the bounded solutions, because one has to eliminate 
the exponential growth of w. Let Q — 3, no — 1. Exponential growth of disappears if 
one of the eigenvalues is zero and two others have opposite signs. For instance, let e 2 = 0, 



ei 



-£3 



1, Ti n) = a n , P ( 



023 



P032 = 0, other Po a p equal i. Then P reads: 
1 i ie-2(*i-M*3) 



ie -(*i+3t2) 1 
j e -2(ti+4t 2 ) ^ 



-(ti+5t 2 ) 



(78) 



which produces the solitary wave solution ii/ -*. 
An open problem remains how to provide physically important reduction w a p = wp a . 
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2.2.2 Pohlmeyer equation 

Eqs. (68,69) corresponding to the eq.(23a) read: 

E tn + E Xn E = 0, (79) 
P tn +P Xn E = 0, n = l,2. 

The general solution of the eqs. (79) is described by the next pair of non-differential equations 
(the appropriate variant of the formulae (70,71)): 

E p = F<f ) {Ep,x 1 -t 1 Ep,x 2 -t 2 Ep), (80) 

Pafi = F^p\Ep, xi — hEp, x 2 — t 2 Ep), 
a = l,...,Q, P = l,.. .,noQ, 

where F^\zq, z\, z 2 ) are related by the system (73), so that only n^Q{Q — 1) of them remain 
arbitrary functions of three arguments. Similarly, eqs. (68,69) corresponding to the eq.(24a) 
read: 

E tl +E Xl E = 0, E t2 +E tl E = 0, (81) 
P h +P Xl E = 0, P t2 +P h E = 0. 

Solution to this system is implicitly described by the next pair of non-differential equations: 

E p = F^\Ep,x 1 -t 1 Ep + t 2 El), (82) 

Pafi = F^p\Ep, xi — tiEp + t 2 Ep), 

a = l,...,Q, 13 = 1,.. .,n Q, 

where F^\zq, z\) are related by the system (73), so that only rioQ(Q — 1) of them remain 
arbitrary functions of two arguments. System (82) corresponds to the particular choice of F^ 
and F( p ) in (80). Namely, 

F^(z ,z 1 ,z 2 ) = F^(z ,z 1 -z z 2 ), F^(z , Zl ,z 2 ) = F^ p \z , Zl - z z 2 ), x 2 = 0. (83) 

Hereafter we consider the "reduced" formulae (80,82) with 

= FW( Zl ,z 2 ), F^(z ,z±,z 2 ) = F^( Zl ,z 2 ), (84) 
FW(z 0:Zl ) = F^( Zl ), F^(z , Zl ) = F^ p \ Zl ). 

On explicit solutions. As it was mentioned above, the explicit solutions correspond to 
Ep = const. Then the matrix P is the only object introducing dependence on x in the field w. 
Let F^p(z) be some smooth function vanishing as \z\ — > oo with single maximum at z — z^ max \ 
The algebraic equations (82) describe two-dimensional surfaces in 3-dimensional space: 

x 1 -t 1 Ep + t 2 E 2 p = z( max \ (3 = l,...,n Q. (85) 

Let F^\z±, z 2 ) be smooth functions vanishing as \z±\ — > oo and \z 2 \ — > oo with single maximum 

at Zi = z( nax \ % = 1,2. Algebraic equations (80) describe two-dimensional surfaces in 4- 
dimensional space: 

xi - hEp = z { ™ ax \ x 2 - t 2 Ep = z^ max \ (3 = 1,..., n Q. (86) 
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Break points manifolds for wave (82). Consider the solutions (82a) with wave profile 
breaking assuming that ti is a physical time and the space variables are (xi,t 2 ). The break 
points manifold Ep is defined by the two requirements. 

The first requirement is that the derivative of Ep in some direction(s) in space (xi,t 2 ) tends 
to infinity. Calculating the partial derivatives 

dzF { p E \z)\ z=Xl _ tlEfj+t2E 2 Ep z F^ E) (z)\ z=Xl _ tlE/3+t2E 2 

Ox^p = ^ ' (j t2 tj p = 2 ' y° ( > 



d= (l + {t l -2Ept 2 )d z F ( f\z)) 



z=xi—tiEp+t2E% 

we see that all of them tend to infinity when the denominator of the above expressions is zero: 

(l + (h - 2Ept 2 )d z F { f\z)) = 0. (88) 

This condition tells us that d = 1 ^ on the surface t\ — 2Ept 2 = 0, so that the only way to 
have the infinite derivative of Ep on this surface is through the infinite derivative of the given 
function F^ E \z) in accordance with (87). To eliminate this case, hereafter we consider only 
those functions F^ E \z) which have smooth behaviour and finite derivatives for z e M. Thus, 
the necessary condition for the break points manifold is 

h - 2Ept 2 ± 0. (89) 

Note that the derivative of Ep is zero in the next direction: 

t-VEp = 0, t= J (E 2 p,-1), V = (d Xl ,d t2 ). (90) 

The second requirement to the break points manifold is that function Ep(xi,ti,t 2 ) must change 
its concavity on this manifold. This requirement is evident in (l+l)-dimensional case and is 
well described in textbooks [18]. It is remarkable, that our (2+l)-dimensional case is very 
similar. To show this, we follow the classical strategy and take Ep as independent variable in 
the neighborhood of the breaking point. There are two possible variants: 

1. x\ is a function of Ep, t 2 and ti, (91) 

2. t 2 is a function of Ep, x\ and t\. (92) 

Consider the case (91). All partial derivatives of X\ are finite in the break points manifold (we 
do not need time- derivative for this analysis): 

\ + {t l -2t 2 Ep)dA E \z) 



d z Ff\z) 



(93) 

z=x 1 -t 1 Ep+t 2 E 2 



d t2 x x = -E 2 a . 

The requirement (88) is equivalent to d E[) Xi = 0. Similarly, the second derivatives of X\ are 
following: 

2t 2 {d z Ff\z)f + dlFf\z) 



(d z Fp \z)) 

d t2 d E xi = -2Ep, 



3 



(94) 

z=xi-tiEp+t 2 El 



>/3 

d 2 t2Xl = 0. 
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We take equation d\ fj x\ = as the second requirement to the break points manifold, which, in 
view of (88), reads: 



(2t 2 - {t 1 -2t 2 E p fd 2 z Ff\z)) 



z=x 1 -t 1 E g +t2E 2 g 



0. 



(95) 



Eqs. (88,95) are supplemented by the eq.(82a). 

It is remarkable, that, starting with eq.(92), we will end up with the same equations (88,95). 
We do not represent details of appropriate calculations. 

The system (82a,88,95) can be solved explicitly: 



dlFf\z) 



(96) 



U=- 



z=Ff\E f) ) 

d z F^\z))\ Ef&FjP(z) 



(E) 



»)■ 



z=F < E \E l3 ) 



x 1 = t 1 E p -t 2 E} + F { p E \E p ), 



where Fp E \z) is inverse of F <yE \z): Fp E ^ ^F^ E \z)j = z. This system of three equations 

describes the break points manifold £p. We can calculate two-component velocity (d^xi, 8^2) 
of the break points manifold differentiating the eqs. (82a, 88, 95) with respect to t\ and assuming 
that variables Ep, x\ and t 2 are functions of t\. 



d tl Ep 



d tl x 1 



2Ept 2 — 1\ 



E p (l2tl + (h - 2E p t 2 fdlF{z)\ z=Xl _ E0tl+E> ) ' 



(97) 



E h 



As a simple example, let 
Then the systems (96,97) yield: 

Ep 



'tlt<2 



F (E \z) = -tanh(z). 



2E, 



to = 



d tl E p 



(El 
(Ej 



l) 2 ' 
l) 3 



3El 



Xi 



Efi-2E\ 



2^(1 + 3£|)' 



duxi 



E/3 

~2'' 



d u t 2 



2E a 



(98) 

arctanh(E /3 ), (99) 
(100) 



Equations (99) describe one-dimensional line in 3-dimensional space. 

The similar analysis of the break points manifolds for the functions P a p shows that these 
manifolds coincide with Ep, i.e. they are described by the same system (88,95) supplemented 
by both eqs. (82). To obtain this result, we need eq.(82b) where P a p, t\ and t 2 are taken as 
independent variable while x\ and Ep are considered as functions of them. 



18 



Break points manifold for wave (80). Consider the wave profile breaking for the solutions 
of the eq.(80a) with time t\ and space variables (x±, x 2 , t 2 ). As usual, the first requirement to 
the break points manifold Ep is that the derivative of Ep tends to infinity in some direction(s) 
in space (x±, x 2 , t 2 ). Let us write all partial derivatives: 



dziFp \z\, Z 2 )\ --i -• i ' i l. , 



d xl Ep — 



d tl Ep 



z 1 =x 1 -t 1 E /3 
Z2=x 2 -t2Ea 



d 



d r „E. 



^z 2 Fp \zi, Z2)\z l =x 1 -t 1 Eg 



Epd Zl Fj 3 E \zi, Z2)\z l =x l -t l Eg 
d 



d t2 Ep 



d 

Epd Z2 F ( j } E \zi, z 2 )\z 1 =x 1 -t 1 E fj 

z 2~ x 2~ t 2 E 



d 



d=(l + hd^Pizu z 2 ) + t 2 d Z2 F { p E \ Zl , z 2 )) 

We see that all derivatives tend to infinity if d = 0: 

[l + t 1 d zl F^ E) (z 1 , z 2 ) + t 2 d Z2 Ff\z x , z 2 )) 



z 1 =x 1 -t 1 Ep 
z 2 =X2-t2Ejj 



zi=xi-tiE 
z2= x 2-t2 E l3 



(102) 



Note that the derivative of Ep is zero in the direction t 



'0, Ep, 1). The second require- 



ment is that function must change concavity in the break points. To relate this requirement 
with the second derivative of some coordinate with respect to Ep we follow the usual strategy. 
Let Ep be independent variable. There are three following cases: 



1. x 2 is a function of Ep, x\, t 2 and t 1 , 

2. x\ is a function of Ep, x 2 , t 2 and ti, 

3. t 2 is a function of Ep, x\, x 2 and t 1 . 



(103) 
(104) 
(105) 



Let us consider the first case in details. One can find the first derivatives of x 2 differentiating 
eq.(80) with respect to Ep, x\ and t 2 (we do not need the derivative with respect to time ti): 



9 E aX 2 



d Xl x 2 = 



1 + hd^Ff'izu z 2 ) + t 2 d Z2 Ff\ Zl , z 2 ) 



d Z2 F { p E) < 



Z\,Z 2 ) 



(106) 



z l =x 1 ~^1 E 3 

z 2= x 2-~ t 2 E {3 



d Zl F {E \ Zl ,z 2 ) 



d Z2 F { p E) i 



Z\,Z 2 ) 



d t2 x 2 = E. 



z 1 =x 1 -t 1 E /3 
z 2~ x 2~ t 2 E [3 



The requirement (102) yelds 8e x 2 = 0. Similar to the previous paragraph, we take the 
condition d\^x 2 = as the second requirement defining the break points manifold. The explicit 
form of this requirement may be found differentiating the eq.(106a) with respect to Ep and using 
eqs.(106) for the first derivatives of x 2 . One gets in result: 



hd Z2 Fp E \zi, z 2 ) ^2(l + hd Zl Fp E \zi, z 2 ))d Zl d Z2 Fp E \zi, z 2 ) — 
t l dlFl E \z 1 ,z 2 )d Z2 F < p E \z l ,z 2 )) - (l + t 1 d Zl Ff\z l ,z 2 )) 2 dlFl E \z l ,z 2 )\ 



(107) 
= 0. 



Z2 =x 2~ t 2 E @ 
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This equation may be given a simpler form using (102) in order to eliminate d Z2 F^ E \zi, z 2 ) (if 
t 2 ^0) ovd Zl F^\z x ,z 2 ) (iftx^O): 

(tid 2 Zl F^\z u z 2 ) + 2t x t 2 d Zl d Z2 F<f\z x , z 2 ) + t\dl 2 Ff\z^ z 2 )) \ = 0. (108) 

z 2= x 2~ t 2 E f3 

The equations (102,108) must be supplemented by the eq. (80a). 

Remark, that starting with the eq.(104) or eq.(105) we will end up with the same two 
conditions defining the break points manifold, eqs. (102,108). 

System (80a, 102, 108) allows one to write the linear system defining the velocity of the break 
points manifold. Assuming that variables Ep, x 1 , x 2 and t 2 depend on t x and differentiating 
eqs. (102, 108) with respect to t x we get: 

(E -x' 1 )d Zl F {E \z 1 ,z 2 ) + (Ept' 2 -x' 2 )d Z2 F {E \z x ,z 2 ) = 0, (109) 
8 Zl F<f \ Zl , z 2 ) - Ept 2 d Zx d Z2 F (E \z x , z 2 ) - E p t 1 d 2 zl F {E \z 1 , z 2 ) + 
x[ (t 2 d Zl d Z2 F (E \ Zl , z 2 ) + t x d 2 Zl F {E \z x , z 2 )) + 

x 2 (t 2 d 2 Z2 F^\z 1 ,z 2 )+t 1 d Zl d Z2 F^\z 1 ,z 2 ^ + 

t' 2 (d Z2 Ff\ Zl ,z 2 ) - Ept 2 d 2 Z2 F {E \z Xl z 2 ) - Ept x d zx d Z2 F {E \z x ,z 2 )) = 0, 

-E f3 t 2 1 ^ Zl F { f\z 1 , z 2 ) + 2t l d 2 Zi Ff\z 1 , z 2 ) - 2E f3 t 1 t 2 d 2 Zl d Z2 F { f\z 1 , z 2 ) + 
2t 2 d zl d Z2 F {E \z x , z 2 ) - E p t 2 2 d Zl d 2 Z2 F i f\z 1 , z 2 ) + 

*i (tld 3 z A E \ Zl , z 2 ) + 2t x t 2 d\d Z2 F {E \z x , z 2 ) + t\d z ,d 2 Z2 Ff \ Zl , z 2 )) + 

x' 2 (tld 2 Zl d Z2 F<f\z u z 2 ) + 2t 1 t 2 d Zl d 2 Z2 F i f\z 1 , z 2 ) + t\d\ 2 Ff \z u z 2 )) + 

4 ( - Ept\d\d Z2 Ff\z x , z 2 ) + 2t 1 d zl d Z2 F^ E \z 1 , z 2 ) - 2E p t 1 t 2 d Zl d 2 Z2 F ( f\z 1 , z 2 ) + 

2t 2 d 2 Z2 Ff \z u z 2 ) - E,t\dl 2 Ff\z u z 2 )) - E'p {^F^(z u z 2 ) + 

Zt\t 2 d 2 zl d Z2 Ff\z x , z 2 ) + 3t x t 2 2 d Zl d 2 2 F {E \z x , z 2 ) + tldl 2 F {E \z x , z*j) = 

where z x — x x — Ept x , z 2 = x 2 — Ept 2 . As a simple example, let 

Ff\z x ,z 2 ) = -tanh(zi + z 2 ). (110) 

Then eqs.(80a,102,108,109) yield: 

x x — —x 2 , t 2 = l — t x , E/3 = 0, (111) 
x\ = -x' 2 , t 2 = -l, E' p = 0. (112) 

Equations (111) describe two-dimensional surface in 4-dimensional space; breaking starts at 
fixed value of Ep: Ep = 0. Similar analysis shows that Ep is break points manifold corresponding 
to the functions P a p (see eq.(80b)) as well. 
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3 New version of the dressing method 



We have described some solutions manifolds to well known S'-integrable models (which were 
called PDEl(ti, t 2 \ w^) in Sec. 2. It is remarkable, that these solutions manifolds cover full 
solutions spaces to the assotiated families of the matrix first-order quasilinear PDEs integrable 
by the method of characteristics (which were called PDE0(t m ; n 0: w), m — 1, 2). 

The purpose of this section is the construction of such dressing algorithm which (a) would 
reproduce the solutions manifold to PDEls available by the classical dressing method and 
(b) would describe the different solution manifold to PDEls in spirit of the method of char- 
acteristics, i.e. we will derive non-differential equations implicitly describing solutions to 
PDEl(ti,i 2 ;«/°)) which are solutions to appropriate PDE0(t m ; n , w), m = 1,2, as well. Our 
consideration is based on the new version of the dressing method, which is subsequent devel- 
opment of the basic ideas of [29, 30, 56, 57, 58]. As a particular result, eqs. (41,61,64) will be 
obtained. 

Let us describe the basic objects appearing in the dressing algorithm. In [30], eq.(3) with 
p = has been derived as a simplest example of the implementation of the dressing method 
based on the following integral homogeneous equation 

J W(A, v\ x)U(u; x)dSl(v) = (A, v\ x) * U(u; x) = 0, (113) 

D 

where \P and U are Q x Q matrix dressing and spectral functions respectively, A and v are 
complex spectral parameters, Q is some measure on the complex plane u, " *" means integration 
over v and x — (x 1 , x 2 , ■ ■ ■ , t±, t 2 ■ ■ ■ ) is a set of independent variables of nonlinear PDEs. 

However, integral equation (113) is not appropriate for our purposes. In order to derive 
PDEls together with assotiated PDEOs, we shell replace the equation (113) with the following 
integral homogeneous equation 

V(\,w,x)*U(v,n;x) = 0, (114) 

where (M + 1)Q x Q matrix spectral function U depends on two spectral parameters; \1/ is 
Q x (M + 1)Q dressing function and kernel of the integral operator. Here integer parameter 
M = dimker Following the strategy of [30], we assume that the solution of (114) is not 
unique but may be represented in the form: 

M 

U(X, fi;x) = J2 U {h '> 3 \\ v\ x) * f®(v, n x), (115) 

3=1 

where f^(v, x), j = 1, . . . , M, are arbitrary Q xQ matrix functions of arguments and 

j = 1, . . . , M, are M nontrivial linearly independent solutions of the homogeneous equation 

(114), i.e 

M 

J2u (h '' 3) {\v,x)*S ij) {v,fi;x) = S {j \v,n-x) = 0, j = 1,...,M. (116) 

3=1 

This assumption causes the single linear relation among any M + 1 independent solutions , 
j = 0,...,M, ofeq.(114): 

M 

C/(°)(A, (i;x) = J2 U (j) (\ v\ x) * F®(v, x), (117) 

3=1 
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where are some Q x Q matrix functions. As we shell see, all are expressed in terms of 
the single solution U through some linear operators either differential or non-differential: 

U U) (X,p;x) = L ij) (X,u) * U(u,fi;x), j = 0, . . . , M. (118) 

Thus, eq.(117) represents a linear equation for the spectral function U. Besides, we will show 
that may be expressed in terms of U using the external Q x (M + 1)Q dressing matrix 
function G(X,p;x), similar to [30]. 

To increase dimensionality of the dressing functions (i.e., the number of variables Xj which 
may be introduced arbitrarily in the dressing functions) and, consequently, to increase the 
dimensionality of the assotiated nonlinear PDEs, we introduce the QxQ matrix function A(X, p) 
and (M + \)Q x (M + 1)Q matrix function A(X, p) by the following generalized commutation 



relation involving 

A(X, v) * ^(z/, p; x) = ^(A, u; x) * A(u, p). (119) 
Similar to eq.(13) we define operators A^ and A* as follows: A^ = A^* ■ ■ ■ * A y , A^ = A* • • • * A . 

j 3 

Consequently, functions A(X, p) and A(X, p) generate the following set of functions 

A {mj) = p {mj) {A), A (mj) = p (mj \A) : p (mj \A)*y = y*p {rnj) {A), (120) 



where scalar functions p( m: >\z) are representable by a positive power series of z, see eq.(9), 
so that A {mj) and A {m ^ are well defined QxQ and (M + 1)Q x (M + l)Q matrix functions 
respectively. 

As usual, x-dependence of the spectral function U appears through the x-dependence of the 
dressing functions. Let ^-dependence of \& be given by the equation 

N 

* tm (A, p; x) + A (mj) (\ u) * * Xj (v, p; x) + C< m > (A, u) *^(u,p;x) = (121) 

9(\,v;x)*C< m \v,VL), 

where A {mj) , A^ do not depend on x; C (m) and are Q x Q and (M + l)Q x (M + 1)Q 
matrix functions respectively, A^ *C- C (m) *A^0, A*C - C * A^ 0. Equations (119) and 
(121) represent the overdetermined linear system for Compatibility condition of (119) and 



(121) yelds: 

pH =A* C (m) - C {m) * A, P {m) = A* C {m) - C {m) * A. (123) 
In addition, equations (121) with different values of m are compatible only if = and 



p( m ) = 0. In other words, if p( m ) ^ and p( m ) ^ 0, then derived nonlinear PDEs will not 
possess commuting flows, at least among PDEOs. 

Now we may obtain the set of different solutions to the homogeneous eq. (114) applying 
A m * and (d tm + Y,f=i A {mj) * d Xj + C (m h) to (114). One gets 

V(X,v;x) *E^ m \v,p;x) = 0, j = 1, 2, (124) 
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where 

E (1 ' m \X,fi;x) = A m {X,u)*U{u,fi;x), (125) 

N 

(X,n;x) = U tm (X,n;x) + J2 ^ (X,u)*U Xi (y,^x) + C<™> (X,v) * U(u,fx;x). 

3=1 

Thus, E^' m ^ are desirable solutions of the integral homogeneous equation. The derivation of 
nonlinear PDEs significantly depends on the value of the parameter M. Hereafter we consider 
M = 1. 



3.1 Simplest degeneration of the ker ^f: M = 1 

In this case eq.(115) reads: 

U(X, n- x) = U^\X, v- x) * fW(v, n- x). (126) 

It follows from the above discussion that we have to introduce so-called external dressing Q x 2Q 
matrix function G(X,/i;x), whose prescription will be explored in Sec. 3. 1.1. Let G be defined 
by the next overdetermined system of linear equations: 

G(X, v\ x) * A(u, fj) = A(X, v) * G(u, //; x) + i?i(A; x)H 2 (n; x), (127) 

N 

G tm (X,fi;x) + J2 G Xj (A, v- x) * A™ (u, fi) - G(X, i/; x) * (u, = (128) 

3=1 

N 

-C (m) {X, v) * G{y, n;x)-J2 T {l) G{X, u; x) * i (mj) (z/, n), 

3=1 

where N is some integer, A and H\ are Q x Q, while H 2 is Q x 2Q matrix functions; functions 
A(mj) _ p( m i)(^4) are defined by the series (27). Require also 

T {m) A = AT (m)^ T (m) Hi = tf^M ( 12 g) 

Functions Hi(X;x) and H 2 (fi;x) are called external dressing functions as well as function 
G(X,fj,;x). The compatibility condition of eqs.(127) and (128) yelds: 

G(X, v- x) * P (m) (v, n) - P {rn) (A, v) * G(v, ft; x) = (130) 

N 

-L^H^X; x)H 2 {n- x)-J2 H i*M> X ) H ^ x ) * A(mj) ^ - H ^ ^)4 m) ^ 2 (/i; x), 

3=1 

where 

P (m) = A * C {m) - C im) * A, (131) 
L^H^X-x) = H lt JX;x) + C^\X,u) * H^x), 

N 

L^H^-x) = H 2t Jn;x) + Y,H 2xj {u;x) * A^\y^) + 

3=1 

N 

J2T U) H 2 {u;x)*A^\u,fi)-H 2 {u;x)*C {m) {u,fi). 
3=1 
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In order to derive the eq.(130), we differentiate the eq.(127) with respect to t m and use eq. 
(128) to eliminate G tm and eq.(127) to simplify the result. After all transformations we end up 
with eq.(130). 

The compatibility condition (130) produces equations defining i = 1,2, in the following 
way. Eq.(130) must be identical in A and p. Solving it, we must take into account that two 
spectral parameters are not separated in the LHS of the eq.(130), while each term in the RHS 
of this equation separates two spectral parameters. We suggest two different solution of the eq. 
(130): 



1. Let 



A*H 1 = 0, (132a) 
pW= p M( i l), V (m) = p {m0) (A), P im) =p (m0) (i), (132b) 
H^ + C^Uh^O, H lx .=0, (132c) 



3 

N N 



H 2tm + H 2 Xj * A {m ^ + TU)H 2 * & mj) - H 2 * = (132d) 

3=1 3=1 

-H 2 *p (m0) {A)*A- 1 , 

where functions p( m °) are defined by the eq.(9). Then eq.(130) becomes: 

G * p (m0) (A) - p (m0) (A)*G = #i# 2 p (m0) (A) * A~\ (133) 

which is equivalent to the eq.(127) in virtue of the eq.(132a). Compatibility of eqs.(128) with 
different m requires p( m ) = and p( m ) = 0. These two conditions provide compatibility of 
eqs.(132) with different m. Otherwise m takes a single value (say, m — 1) and assotiated 
nonlinear PDEs have no commuting flows, at least among PDEOs. 

Remark. Eq.(133) requires that positive series for p*™ ^) start with z 1 , which corresponds 
to a( m0; °) = in eq.(9). However, it will be shown in the Sec. (3. 1.2) (see a paragraph after the 
eq.(147)), that this is not a strong restriction. 

2. Let 

A * #i ^ 0, C (m) = C im) = 0, C (m) = 0, (134a) 
H ltm = 0, H lx .=0, j = l,...,N, (134b) 

N N 

H 2tm + H 2 Xj * & (mi) + E tU)h * * A{mj) = °- ( 134c ) 

3=1 3=1 

Appropriate nonlinear PDEs possess commuting flows since the eqs. (128) with different m are 
compatible, as well as eqs. (134) with different m. 

It is evident that the eq.(122) is satisfied for both cases (132) and (134). 

3.1.1 Spectral system for U(X,fi;x) 

Now we are ready to derive the linear spectral system. Eq. (126) means that any two solutions 
of the homogeneous equation (114) are linearly dependent. In particular, any solution E^' m ^ 
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(125) is linearly dependent on U, i.e 

£ {1;1) (A,/i;x) = U(X,is) *F(v,n;x), => (135) 
A(X, v;x)*U (v, (j,; x) — U (A, v\ x) * F(v, (/,; x), 

E (2 ' m) (X,fx;x) = U(X,is) *F (m) (z/,/i;x), => (136) 

N 

U tm (X,fi;x) + J2 A (m3) (A, v) * U Xj (u, n x) + ( A, u) * U(u, x) = 

3=1 

U(X,is;x) * F {m) (u,fi;x), m = l,2,.... 

Eqs. (135,136) represent the overdetermined linear system (the spectral system) for the spectral 
function U(X,fj,;x). Remember that solution of eq.(114) is not unique. To obtain uniqueness 
we introduce one more equation for the spectral function U. This equation can be largely 
arbitrary, but, in order to derive the simplest nonlinear PDEs, we select the following equation: 

G{\,v;x)*U(y,n;x) = I8{\- n), (137) 

so that U is the unique solution of the system (114,137). In other words, the equation (137) 
fixes the function /^(A, x) in the eq.(126). Now, applying G* to the eqs. (135, 136), one gets 
the following expressions for F and F^\ 

F(is,(i;x) = G(\,u;x)*E {1 ' l) (u,fi;x) = (138) 
A(A,/i) + H 1 (\;x)H 2 (v;x) * U(u : /j,;x) : 
F (m) {y, At; x) = G(X, v\ x) * E {2 ' m) {y, n\ x) = 

N 

(7 (m) (A,/i) + ^T (J) G(A, v\ x) * (u, v) * U{u,ji]x) + 

3=1 

N 



(g(\,v;x) * A {mj \u,u) *U(v,n;xf) , m=l,2, 



3=1 

Thus, eqs. (135, 136) become the overdetermined spectral system for the spectral function U (A, x). 
We consider the particular examples of the spectral systems and of the appropriate nonlinear 
PDEs in Secs.3.1.2 and 3.1.3. 



Spectral functions depending on single spectral parameter; fields and remarkable 
reductions. Eqs. (135,136) depend on two spectral parameters due to the spectral function 
U(X,/i;x). However, functions of single spectral parameter appear in the algorithm naturally. 
These functions are following: 

V U) (X;x) = U(\,n;x) *A j (jj, v) *H 1 (v;x), (139) 
WV\n;x) = H 2 (X;x)*A j (X,u)*U(u,i2;x). 

They satisfy the spectral equations with single spectral parameter which appear after applying 
*A* * H-l and H 2 * AH to the eqs.(135,136), see Secs.3. 1.2,3. 1.3. 

The dependent variables (or fields) of the nonlinear PDEs are expressed in terms of the 
spectral and dressing functions by the next formulae: 

w {kn \x) = H 2 (X; x) * A k (X, /i) * V {n \fi; x) = W {k \fi; x) * i n (/i, v) * H x {y- x). (140) 
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This definition of the fields suggests us the following two types of reductions: 

no - 1 

1. A n ° * H 1 = & * H ir {j) (141) 

3=1 
no — 1 

w (kn ) = J- W (^\ VJfc, k= °4 32) (35), 

3=1 
fco-1 

2. # 2 * A ka = Y r U) H 2 * A j (142) 

3=1 

w (k n) = J- r Ci) w C*0 Vn n= ^ 34) (36). 

3=1 

It will be shown in Sec. (3.2.2) (see text after the eq.(205)), that must be scalar constants 
(remember that these parameters are matrices in the algebraic approach, see the paragraph 
after the eqs.(21)). We consider only the reduction (141). 

In accordance with Introduction, we refer to the nonlinear PDEs corresponding to the 
reduction (141) as PDEOs, while the nonlinear PDEs without any reduction will be referred 
to as PDEls. Although the nonlinear PDEs may be derived for any given functions p^ m ^ and 
p(mj)^ ^ j g difficult to write the nonlinear PDEs keeping these functions unfixed, in general. 
The example considered in Sec. 3. 1.2 is an exception. 

3.1.2 Case (132). Quasilinear first order matrix equations solvable by the method 
of characteristics. 

In this section, we briefly reproduce results obtained in [29, 30] by a different method. First 
of all, remark that the eq.(132a) is nothing but the reduction (141) with no = 1. Thus, the 
derived equations will be PDEOs. In this case eqs. (135,136) yield: 

A*U = U*A + U* E X E 2 * U, (143) 

N 



u tm + Y A(mj) * u Xj + c {m) *u = u* c (m) + 

3=1 

N N 

T ij) (A (mj) + H X H 2 * A (mi) * A' 1 * f/) + ^(#i#2 * A imj) * A' 1 * U) x ^j 



3=1 3=1 

Applying *Hi to the eqs. (143) one gets (V = V^°\ w = -u/ 00 )) 

A*V = Vw A imj) * V = Vp (mj) {w), (144) 

N 

VtjX; x)+J2 A {mj) (X, v) * V Xj (v; x) + C (m \X, u) * V(v; x) = (145) 



3=1 

N N 

V(X; x) T {j) P {mj) (w) + V(X; x) ^ p^\w). 

3=1 3=1 
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or, substituting the eq. (144) into the eq.(145), one has: 



N 



VtJ A; x) + J2 V Xj (A; x)p™ (w) + C< m > (A, v) * V{y; x) = (146) 

3=1 

N 

V{X;x)Y,T U) p {mj) {w). 

3=1 

After applying H 2 * to this equation one ends up with PDEO: 

N N 

w tm + w *>P {mj) H + P (m0) H = \™> Yl TU) P imj) HI • (147) 

3=1 3=1 

It was mentioned above, that this equation has commuting flows only if p( m0 ) = 0. Then 
m = 1, 2, . . . in this equation. Otherwise, m must be fixed, say, m — 1. 

In accordance with the Remark given after the eq.(133), power series for p^\z) must start 
with z 1 . Otherwise, multiplying eq.(147) by w from the right and replacing t\ — > xat+i we 
result in the PDEO having the form of "stationary" eq.(147) satisfying the condition of the 
above remark. 

Eqs. (144) and (145) has to be taken as the spectral system with the spectral function 
V(X; x) corresponding to the nonlinear PDE (147). In [29] this equation has been derived with 
fixed m = 1, p^ = p^\ j > 0, p^ = p, p^ = pb\ 3 . In addition, p = in [30], where another 
variant of the dressing method has been introduced. However, both techniques developed in 
[29] and [30] may be applied to the eq.(147) as well. 



3.1.3 Case (134). PDEls and assotiated PDEOs: iV-wave and Pohlmeyer equa- 
tions 

The main feature of the eqs. (134) is that they set = and do not force the reduction 
(141). Nonlinear PDEs corresponding to this case give rise to PDEls. Then, been imposed, 
reduction (141) reduces PDEls to PDEOs. 

We study nonlinear PDEs corresponding to the particular choice of the functions p( mj ) and 
p(mj) m For instance, let = s im) A5 mj , A (m ^ = A5 mj , s (m) are scalar constants. Then the 

spectral system with two spectral parameters, eqs. (135,136), reads: 

A*U = U * A + V {0) W {0) , (148) 
U t +s (m) A*U x = U*AT (m) + V i0) T (m) W (0) + s (m) V (0) Wi°\ m=l,2. 

The appropriate spectral system for the spectral functions V^(X;x) with single spectral pa- 
rameter appears after applying *A n * Hi, n = 0, 1, . . . , to the eqs. (148): 

A(X, v; x) * (u;x) — y^ n+1 ^(A; x) + V^°\X; x)w^- 0n \x), (149) 

v£>(A; x) + (a(\, v) * V^{y; x) - V^{X; x)w^\x)) = V^ n+1 \X; x)T™ + 

V (0 \X;x)T {m) w {0n) (x), m = 1,2. 

In principle, applying H 2 * A k * to the eq.(148), one could get the linear equations for another 
set of functions of single spectral parameter W^ k \p; x). However, these equations will not be 
used, so that we do not represent them here. 
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The system of nonlinear PDEs may be derived applying H 2 * A k * to the eq.(149). 

w ((k+l)n) = w (k(n+l)) + w (k0) w (0n)^ (150) 
w [ k ") + S M( w ((fe+l)n) _ ^(^^(On)) = (151) 
w (k(n+l)) T (m) _ T (m) w ((k+l)n) + w (kO) T (m) ^(On) ^ TO = ]_ 2. 

Using eq.(150) to eliminate u;(( fe+1 ) n ) from the eq.(151) and putting /c = in the result we get 
(see eq.(32) for definition of w^): 

wt2 + s (m) (^fc +1) + w ( xlw in) ) ~ [w (n+1 \ T (m) ] + [T {m \w {0) }w in) = 0, (152) 
n — 0,1, ... , m = 1,2. 

Fixing n = and eliminating from this system of two PDEs, we get PDEl^!, t 2 ; w^). 
Reduction (141) yields PDE0(t m ; n , w), m — 1,2: 

w tm + s^w Xm w+[Tif,w}w = 0, m = l,2. (153) 

From another point of view, eliminating ty( fc ( n+1 )) from the system (151) and putting n = 
one gets: 

+ s M (*t +1) - w^w^l) + [T^ m \w^] + wW[w<®,tW] = 0, (154) 
A; = 0,1,..., m=l,2, 

where fields are defined by the eq.(34). Fixing k = and eliminating from this system 
of two PDEs, we get the same PDEl(ti, t 2 ; w (0) ). Reduction (142) yields PDE0(t m ; k ,w), 
m = 1,2: 

w tm + s™ww Xm + w[w, T fc ( o m) ] =0, m = 1, 2. (155) 

iV-wave and Pohlmeyer equations with apporpriate PDEOs and spectral problems follow from 
eqs.(149-155) with s (m ^ = (see Sec.2.1.1) and = 0, s (m ^ = 1 (see Sec.2.1.2) respectively. 

Eq.(154) can be formally obtained by the transposition of the eq.(152) with replacements 
w T — ► — w, — > — T^ m \ similarly to the relations between eqs.(33) and (31) and between 
eqs.(58) and (57). Because of this similarity, we will study only eqs. (152,153). 

3.2 Solutions space described by the dressing method 

The dressing algorithm describing solutions space to PDEls consists of two basic parts: 

1. The non- evolutionary part, where we do not need the particular dependence of the dressing 
functions on X{ and tj. As a result, we derive the system of (non-differential) equations 
which is valid for any possible ^-dependence of the dressing functions. 

(a) Solve the eq.(119) as an equation defining 

(b) Solve the spectral equation (148a) as an equation establishing the restrictions for 
the spectral function U. 

(c) Solve the eq.(127) as an equation for G. 

(d) Solve the eqs. (114) and (137) as equations for U taking into account result of n.lb. 
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2. The evolutionary part, where we presume the special x-dependence of the dressing func- 
tions. 



(a) Solve eq. (121) as an equation defining x-dependence of the internal dressing function 

(b) Solve one of the systems (132) or (134) defining Hi, i = 1,2. This step defines 
x-dependence of the external dressing function G(X, /x; x) owing to the eq.(127). Eq. 
(128) will be automatically satisfied. 

Details of the dressing algorithm producing the classical solutions manifold to PDEls are given 
in Sec. 3. 2.1. In order to describe solutions manifold to PDEls assotiated with reduction (141) 
one has to evaluate the same algorithm with eqs.(114, 137,148a) modified by the reduction 
(141), see Sec. 3. 2. 2 for details. 

Before describing all steps of the algorithm in details, we give a description of the particular 
representation of the dressing and spectral functions which is necessary hereafter. First of all, 
in order to fit the requirement M — 1, we introduce the next block-matrix representation of 
the functions: 



V(\,[i;x) = [tp (X, fi; x) Vi(A,//;x)], 



u (X,/jl;x) 
ui(\,fj,;x) 



V {j \X-x) 



,(i) 
Si) 



(X;x) 
(A;x) 



U(\,fi;x) = 



(156) 



G(X, n; x) = [g (X, /x; x) gi(X, /x; x)\, 
Hi(X; x) = h±(X] x). 



H 2 (X;x) = [h 20 (X;x) h 2 i(X;x)], 



Any function in the RHS of formulae (156) is Q x Q matrix function, so that G and H 2 are 
Q x 2Q matrix functions; U and are 2Q x Q matrix functions; and Hi are Q x Q 

matrix functions. We also have to fix the functions A(X, fj), A(X, /x) and A(X, /x) (compare with 
eq.(52)): 



^(A, fj) = A(X, fj) = X5(X - /x)/, A(X, /x) = X5(X - /x)/ 2 . 



(157) 



3.2.1 Dressing algorithm: classical solutions manifolds to PDEls 

Non-evolutionary part of the dressing algorithm. We describe all steps of the non- 
evolutionary part of the algorithm given in the beginning of the section (3.2) for the case (134) 
without reduction (141). 

(la) Function \l> is defined by the eq.(119), which reads 

(A - /x)*(A,/x;x) = 0. (158) 

Its solution is following: 

*(\,H;x) = 4>(\;x)6(\ - fi), *(A; x) = [ip (X; x) MX;x)], (159) 

where fy, i = 0,1, are Q x Q matrix functions. 

(lb) As it was mentioned above, the main feature of the spectral system described in this 
paper is the presence of the spectral equation which has no derivatives with respect to x, see 
eq.(148a): 

U(X,fi;x)(X-fi) = V m {X;x)W m {fi;x). (160) 
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This suggests us the next representation for U : 

(161) 



y(A,|x;x) = V{0) ^ )W{0) ^ ;X) +U (X;x)6(\-^ U (X;x) 

A — jl 



u 00 (X;x) 
u 01 (X;x) 



where uoi, i — 0,1, are Q x Q matrix functions, 
(lc) Similarly, eq.(127) reads 



G(X,fx;x)(X- fx) = -H 1 (X;x)H 2 {n;x) (162) 
which suggest us the next formula for G: 

G(X, n x) = _ H ^*) H *<<K g) + G (X; x)S(X - jx), G (X;x) = [g 00 (X;x) g i(X;x)}, (163) 

A — jl 

where g 0i , i — 0,1, are QxQ matrix functions and G (A; x)5(X — /i) is a solution of the eq.(128). 
(Id) Next, substituting the eqs.(161) and (163) into the eq.(137) one gets 

nn x /" , H 2 (u;x)V^(u;x)W^(fi;x) ^ H 2 (fi; x)U (fi; x) 

-ifi(A;ar) / di/ r ffi A;x + 164 

J (A - v){y - /i) X-pL 

D 

G (X;x)V^(X;x)W^(/i;x) „ STT .... , TS/ . 
X — ij, 

This equation must be identity for any A and ji, which suggests us to split the eq.(164) into 
two matrix equations. The first equation reads: 

G (X;x)U (X;x) = I, (165) 

which is the system of Q 2 scalar equations for 2Q 2 elements of Uq, i.e. underdetermined system 
of scalar equations. The second equation reads: 

1 [E 1 {X;x)W{n;x) + H 1 {X;x)E 2 {n;x)] = 0, (166) 



A — ji 

E X {X- x) = G {X; x)V^{X; x) - H^X; x) ! dv H ^ x ) v{0) ^ x ) _ H ^ x . ^ 

J X — v 

D 

x f r H 2 (v,x)V^(v,x)W^U;x) TT , , TT , , , 

D 

In turn, the eq.(166) must be splited into two matrix equations: 

E 1 (X;x) = 0, (167) 
E 2 (/i;x)=0. (168) 

Due to our choice of the last term in the expression for E 2 , eq.(168) coincides with eq.(161) 
after application H 2 * to it, which is necessary condition. The last term in expression for E 1 
compensates the last term of E 2 after substitution in eq.(166). 

The matrix equation (167) represents Q 2 scalar equations for 2Q 2 elements of the matrix 
function V^°\ while the matrix equation (168) represents Q 2 scalar equations for Q 2 elements 



30 



of the matrix function . Thus, matrix equation (167) is underdetermined systems of scalar 
equations. 

The rest of equations for elements of and Uq follows from the eq.(114) after substitution 
the eq.(159) for \l/ and the eq.(161) for U: 

(A; x) y( ° )(A; f W(0) ^ x ^ + i( A; x )Uo(X; x)5(X - = 0. (169) 

Since eq.(169) must be identity for any A and //, it must be splitted into the following equations 
for and U : 

^(X;x)V {0) (X;x) = 0, (170) 
tf(A;a;)E/o(A;a;) = (171) 

Each of these matrix equations represents Q 2 scalar equations for 2Q 2 elements of the matrix 
functions and Uq respectively. Thus, the system (165,167,168,170,171) is the complete 
system for elements of V^(X; x), W^(fj,; x) and Uq{X\ x). Having these functions, the function 
U(X,/i;x) may be constructed using the formula (161). This ends the non-evolutionary part of 
the dressing algorithm. 

All in all, in order to construct the spectral function U, one has to solve 

1. the eqs. (165) and (171) for U , (172) 

2. the eqs. (167) and (170) for V^°\ 

3. the eq.(168) for W {0 \ 

After the function U has been constructed, one must use the definition (140) in order to 
construct the solution to PDEl(ti, t 2 ] w^), i.e. the function k/ 00 ) = w^°\ However, it is simple 
to observe, that essentially important for construction of — H 2 * are eqs. (167) and 
(170) defining V^. Moreover, eq.(167) is equivalent to the classical <9-problem for Self-dual type 
S'-integrable equations. Let us write this equation in standard form. To simplify derivation, we 
take 

H 1 (X;x) = I, (173) 
which is consistent with eq.(134b). Eq.(170) yelds 

v^^x) = -^\X;x)MX;x)v^(X;x). (174) 
Then eq.(167) gets the next form: 

x)4 0, (A; *) + f iv ^lM^l = /, (175) 
J X — v 

D 

where 

0(A; x) = g 00 (X; x) - g m (X; x)4>^ l (X; x)ip (X; x), (176) 
X(X; x) = -h 20 (X; x) + h 21 (X; x)^r" 1 (A; x)4) (X; x). 
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Introduce the new spectral function 

v(X;x) = <f)(\;x)v { o\\;x). (177) 



Then eq.(175) yelds 



J dvR(X, v\ x)v{v- x) = I, (178) 



D 



where 



R(X, v\ x) = ^ ^ + I5(X - u), R(v;x) = x(v,x)(f)-\v,x). (179) 
A — v 

Here R{y\ x) is a new dressing function and R is a new kernel of the integral operator. By 
construction, dimkeri? = 0, so that eq.(178) is uniquely solvable for v. Linear PDE for R 
follows from the linear PDEs for the functions ipi, goi and h 2 i, i — 0,1. These PDEs are 
eqs.(121), (128) and (134c), which read (taking into account eq.(134a)): 



(&(A;x)) +^p^)(A)(^(A;x)) =0, (180) 

N N 

(g 0i (X;x)) +X;P (wu " } (A)(flbi(A;x)) + ^T^)(A) g 0i (X;x) = 0, 
m j=i j j=\ 

N N 

h 2i (X;x)) +J2p (mj) (X) (h 2i (X;x)) +J2T U) p {mj) (X)h 2l (X;x) = 0, (181) 

' . * ' Xn 



J' = l 

TV 

AT AT 

(fc2i(A;x)) + £ 
where i = 0, 1. From this system, we obtain the linear PDE for R: 

N N 

R t jX;x) + J2p {mj) WR* j ( X -> x ) + Y,P {mj) W [T®,R(\;x)]=0, (182) 

3=1 3=1 

Eqs.(178,179,182) represent the classical <9-problem for PDEls, [52, 53]. Eq.(140) for field w (0) 
reduces to the next one: 



w (o) 



(x) = J dX R{X;x)v(X;x). (183) 



D 



Evolutionary part of the dressing algorithm. In this case, nn.2a,2b of the dressing 
algorithm reduce to the solution of the single eq.(182): 

R(X;x)= (184) 

N ( 1 \ N 2 AT 2 

<ig e J=1 v m=1 y J'= lm=1 -Ro(A, g)e J=lm=1 



where i?o is arbitrary Q x Q matrix function, q = (g 1? . . . , g^), parameters g« are complex in 
general and T> is some integration region in space of vector parameter q. 
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iV-wave and Pohlmeyer equations correspond to p^ m ^ = 0, p( mj \\; x) = XS m j and p^ m ^ = 0, 
p( mj \\;x) = \5 m j respectively. 

The evident manifold of explicit particular solutions corresponds to the following choice of 

R : 

(Ro(X;q)) a p = < j , a,p = l,...,Q 

{ 0, a = 

with scalar parameters aa\a 2 ^ o!ila A V i,j and (0:1,0:2) 7^ (03,04) and arbitrary functions 
c^l(q). Such solutions have the form of fractional rational function of the exponential func- 
tions with linear in x arguments. We do not consider these solutions in details. Instead, we 
concentrate on another manifold of particular solutions in the next subsection. 



3.2.2 Dressing algorithm: reduction (141) and appropriate solutions to PDE0(t m ; n , 
m = 1, 2, and PDEl(ti, t 2 ; w^) 

Non-evolutionary part of the dressing algorithm. We represent the non-evolutionary 
part of the dressing algorithm describing the solutions manifold to both PDE0(t m ; no, w), m — 
1,2, and PDEl(t 1? t 2 ; w^) assotiated with reduction (141). Solution of (141a) should be taken 
in the next form: 

no 

H 1 (^x)=Y t h^(x)6( l jL-b j ), (185) 
3=1 

where h^\x) are diagonal matrix functions and bj, j = 1, . . . ,n are scalar constants solving 
the next linear algebraic system: 

= ^b)r®, (186) 
i=i 

where the parameters are scalars. Requirement to have scalar 6j and, as a consequence, 
scalar is related with transformation from the eq.(204) to the eq.(207), see text after the 
eq.(205). Eq.(186) has n different roots bj, j = 1, . . . ,n , which agrees with the summation 
limits in eq.(185). 

We describe the manifold of the particular solutions to PDEl(ii, t 2 ; w^) corresponding to 
the reduction (141) starting with the original system (114,137,148a) modified by the reduction 
(141) as follows. 

Applying *A j * Hi to (114), j — 0, . . . , n — 1, one gets: 

^*V = 0, (187) 

where V is 2Q x n^Q matrix function (47). Similarly, applying *A : > * Hi to the eq.(137) one 
gets 

G(\,v;x)*V(i>;x) = P(\;x), (188) 

where P is a row of n blocks 

P=[Hi A* Hi ■■■ i" 0_1 *H 1 \. (189) 



33 



Finally, the spectral equation (148a) can be written in the next form after applying *A j * Hi. 

A{\v)*V {] \v-x) = V u+1) (X;x) + V {0) (X;x)w^\ (190) 

or 

A(\,fj,)*V(fji;x) = V(\;x)w(x), (191) 

where w is UqQ x n^Q block matrix (21a). 

In order to solve the system (187,188,191) we assume diagonalizability of w, i.e. w is 
representable in the form (50). The case when w is not diagonalizable but reducible to Jordan 
form is more complicated ind will be discussed in different paper. Now we describe all steps of 
the non-evolutionary part of the dressing algorithm using A, A and A given by the eqs.(157). 

(la) coincides with the previous case, see eqs.(158,159). 

(lb) Eq. (191) after multiplication by P from the right yelds the eq.(53), which we write in 
the next form: 



noQ 



^2y a7 (fi;x)P l(3 (x) =V a p{x)5(n-Ep), 

7=1 

a = 1, ... ,2Q, /? = !,.. . ,n Q, 



(192) 



where V has the structure given by the eqs.(54) and 



V^{x) 



•o 



Vq(X) 



v\ (x) 



, % = 0, . . . ,n - 1. 



(193) 



(lc) coincides with the previous case, see eqs.(162,163). 

(Id) Multiplying the eq.(187) by P from the right, using the eq.(159) for \1> and the eq.(192) 
for V one gets: 



2Q 1 2Q 

J^^ ai (E p ;x)V lP (x) =^Y t $ j )c«(Ef,x)(y j )^(x) = 0, 

7=1 j=0 7=1 

a = l,...,Q, (3 = 1,.. .,n Q, 



(194) 



where 



V, = [v 



(o) 



V 



, 3 = 0,1. 



Eq.(188) reads: 

G(X, v- x) * V(i/; x) = Hi(X; x)P(X), P(X) = [I XI ■ ■ ■ A n ° -1 /]. 



(195) 



(196) 



Multiplying this equation by P from the right, substituting G from the eq.(163) and V from 
the eq.(192) one gets: 



2Q 



{H 1 ) a (X;x)[-J2 



7=1 



(H 2 ) a ^(E ;x)V 1 f 3 (x) 

X — Ea 



J2K(X)P l( s(x) + 



(197) 



7=1 



^2g 0j (X;x)(V j (x)) afS 5 (X - Ep)) = 0, a = l,...,Q, (5 = 1, . . . , n Q. 

3=0 
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Eq.(197) consists of terms with separated spectral parameters. Thus, it must be splitted into 
two equations: 

mafrx) xl Ep =0, (198) 

1 

^2g Oj (E p ;x)(y j (x)) a p = 0, a = l,...,Q, (3 = 1, . . . , n Q. (199) 

3=0 

Consider Hi defined by the eq.(185). Then the equation (198) becomes equivalent to the next 
one: 

1 2Q n Q Q 

E E faUiEfi x ) ( E P ~ hY 1 (YiMx) - E Pcn(bi)Pyp(x) = 0. (200) 

3=0 7=1 7=1 

Eqs. (194, 200) must be viewed as a system of equations for V, while eq. (199) is a matrix 
equation for elements of E and P. These three equations may be written in terms of w in the fol- 
lowing way. First of all, the eq.(194) relates V,-, j = 0, 1: Vi(x) = -^ 1 {E^x)^ {E p ;x)Vq{x). 
Substituting this relation into eq.(199) and into eq.(200) multiplied by (Ep — bj) from the right 
one results in 

2Q 

J^(^) aj (E p ;x)(V (x)U = 0, (201) 

7=1 

a = l,...,Q, f3 = 1, . . . ,n Q. 

Q noQ 

J2h^(Ep;x)(V ),p(x) = E^ T (^)^W(^ - bj), (202) 

7=1 7=1 

a = l,...Q, f3 = 1, . . . ,n Q, j = l,...,n , 

where 

^(-)(A; x) = g 00 (X; x) - g 01 (X; x)^ 1 " 1 (A; x)^ (A; x), (203) 
h(~\\;x) = h 20 (\;x) - h 2 i(\; x)^ 1 (\; x)i/> (\; x) . 

Next, eliminating V from the eq. (201) using the eq.(202) one gets 

Q noQ 

E E R ^( E ^ x)Ep 1 P, 112 (b j )P 12 p(x)(Ep - bj) = 0, (204) 

7i=l 72=1 

a = l,...Q, (3 = 1,.. .,n Q. 

where 

R{\; x) = ^ (_) (A; x){h { - ) )- 1 {\; x)X. (205) 
Finally, applying P -1 from the right, using the fact that bj are scalars and the evident identity 

P(bj) (w - bjl) = P(0) w (206) 
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one gets the resulting equation in terms of w: 
Q 



7=1 



R aj (w;x) = 0, a = \,...Q, (3 — 1, . . . , n Q, (207) 

J 7/3 



which is an equation for the elements of Q x Q matrices w^ 0j \ j — 0, . . . , n — 1. This ends the 
non-evolutionary part of the dressing algorithm. 

We see that, along with the solution U(X,fj,;x) to the spectral problem, we have described 
implicitly the function w(x) (solution to the apporpriate PDEO) and, as a consequence, the 
function w^°\x) (solution to PDE1). As a simple example, let no = 1. Then A * Hi = 0, 
P(k) = I,w = w m . The eq. (207) gets the next form: 



7=1 



Q 

= 0, a,(3 = l,...,Q. (208) 

J 7/3 



R a7 (w; x) 



Eq.(207) holds for any PDE1 and assotiated PDEO. In order to show the equivalence of 
this equation to the non-differential equations derived in Sec. 2 (see eqs.(41,61,64)), we must 
introduce the particular dependence of R on variables and x iy which will be done in the rest 
of Sec.3.2.2. 

Evolutionary part of the dressing algorithm. As we have seen, the structure of eq.(207) 
does not depend on the particular x-dependence of the functions if)^ and defined through 
the x-dependence of the functions G and H 2 (see eqs. (121,128,132,134)). Here we describe 
the ^-dependence of (207), related with the x-dependence of \I/ (eq.(121)), G (eq.(128)) and H 2 
(eqs.(132d) or (134c)). Remark, that evolution of Hi given by the eqs. (132c) or (134b), is not 
needed in resulting formula (207). Emphasise that index m appearing in the equations of this 
section is meaningful only in the case (134), when V {m) = P (m) = and = 0. Otherwise 
this index must be fixed indicating that the appropriate PDEs do not possess the commuting 
flows. 

Functions and are defined by the eqs.(123,131,132b): 

C^(X,im) = C7( m )(A,u) = -(p(A)5(A-u)) A J, (209) 
C( m )(A,u) = -(p(A)5(A-u)) A / 2 . 

Eqs.(121) (in view of eq.(159)) and (128) yield the following equations for the functions ipj and 
g 0j respectively (j = 0, 1): 

N 

(A(A;x)) 4 +Y,P (mj) W(M^x)) x -(^(A;x)p( m °)(A)) A = 0, i = 0,l, (210) 

j=i X] 

N 

(g 0i (\-x))^+J2p {mj) W(9oi(X-,x)) x - (^(A;x)p( m °)(A)) A + (211) 
j=i X] 

N 

J2T^p^\X)g 0l (X;x) = 0, i = Q,l, 

3=1 
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Combining these equations and taking into account the definition of ^ ) ( eq.(203a)) we obtain 
equation for ip^: 

4Z ] (A; x) + P {m3) W1>£ (A; x ) - (V^ (A; x)p^ (A)) ^ + (212) 
i=i 

TV 

^ T 0')p("y)( A ) ^(-)(A; X ) = 0, % = 0,1. 
j'=i 

Mostly general equation for if 2 is eq.(132d), which may be written in the next form in view 
of eq.(132b): 

N N 

H 2t m + H ^ * A(mj) + Yl T(j)H 2 * A (mj) ~H 2 * A' 1 * C {m) * A = 0. (213) 

3=1 3=1 

Then, in particular, eq.(134c) corresponds to = 0. Introduce the function H(X;x) by the 
formula: 

H 2 =H 2 * A. (214) 
Then the above equation takes the next form after applying *A~ l : 

N N 

K2 tm + n ^ * A(m3) + Yl t0) ^ 2 * A(mj) - n * * cim) = °- ( 215 ) 

3=1 3=1 

Function H 2 has the structure analogous to H 2 : 

H 2 (p; x) = [h 20 (p; x) h 21 (p;x)}. (216) 

Then eqs.(215) yields: 

N 



(h«(A;x)) +J>^)(A) (h«(A;x)) + (217) 

3=1 

N 

J2 Tij) P imj) ( X ) h 2,(A;x) - (h 2l (A;x)) A p( m °)(A) = 0, i = 0,1. 



Combining eqs.(210) and (217), taking into account definition of (203b) one gets 

N 



h[-\\;x) +Y,P {mi) (\) h^ XJ (X;x) + (218) 

3=1 

N 

T^pW (A) h<-> (A; x) - (h^ (A; a;)) A P (m0) (A) = 0, 2 = 0,1. 



j'=i 
where 

h ( " } (A;x) = h 2O (X;x)-h 21 (X;x)^ 1 (X;x)^ (X;x), h^\\; x) = h (_) (A; x)X. (219) 
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Finally, linear PDE for R(X;x) defined by the eq.(205) reads: 



(220) 



N 

J2~P (mj) (\) [T^,R(X;x)]-^R(X;x)p^°\X)) x = 0, i = 0, 1. 
i=i 

Thus, evolutionary part of the dressing algorithm, nn.2a,2b, is reduced to the solution of the 
single eq.(220). We describe two particular examples. 

Solutions to the PDEOs (147). In this case no = 1 so that w = w^°\ It has been shown 
in Sec. 3. 1.2, that index m must be fixed. Let m — 1. Use the following notations: 

t = t u p U) = p {ll \ p U) = p (l3) . (221) 

In general, when p^ ^ 0, the eq.(220) may be integrated by the method of characteristics. 
Eq.(220) yelds: 

- = -p (0) (A), 



dt 



(222) 



A? 



Integrating the eqs.(222) in quadratures we receive: 

A 

t = - 



dX 



1 



pW(A) 
dt 



+ Cj, j = 1, 



(223) 



,N, 



R = e 



l IV l iV 

■fdtE TV)pU)(\(t))C(co, ...,C N ) fdt E TO)pO)(A(t)) 

P (0) (A) 6 



where q are integration constants and C(z , • • • , zjv) is an arbitrary function of iV+1 arguments. 

Integration constants q = q(x, x), i = 0, . . . , N, must be expressed in terms of A and x 
using (223a,b). Let 

C(z Q , ...,z N )=zo- F{z , ...,z N ). (224) 

Then the eq.(207) yields 

(c (w,x)) aP = (225) 



Q 



7=1 



2 = 1 



^F a7 (c (tu; x), . . . , cat(u>; x))e 



/*Er 7 0) pO)(A(t)) 

2 = 1 



7/3 



a = 1, . . .Q, (3 = 1,.. .,n Q. 

This equation describes the solutions space to the PDE0 (147) with m = 1 and notations (221). 
Of course, the same result may be obtained using the algebraic approach. 
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Construction of the solutions manifold to the PDEl(ti, t 2 ; w^) and PDEO(i m ;n ,w), 
m = 1,2, derived in Sec.3.1.3. If p (m0) = (see eq.(134)) then eq.(220) may be solved 
explicitly: 

JV / 2 \ JV 2 JV 2 

R(X;x)= dqe^ 1 ^ m=1 y *= lm=1 R (\, q)&= lm=1 ,(226) 

where q — (q±, . . . , gAr), parameters are complex in general and integration is over some region 
of the space of parameter q. The proper choice of i? is following: 

R (\,q) = \-F(\,q). (227) 

Then the eq.(207) gives us 

^r-E E Ti J V m3) (to)** 



2 JV 

r - E E Ji"c v " j 'wtm / x , 

w a/3 = 2Z[e Far^an/no-JJp^H^,..., (228) 



7=1 m=l 

2 JV 



m=l 

a = 1, . . /3 = 1, . . . ,noQ, 



J 7/3' 



where we use notation 

JV 

F(zq,...,z n ) = J dqF(z , qi ,...,q N )e^ qjZj . (229) 

Solutions to iV-wave equation. If p( m ^ = 0, p^ m ^ = A5 m j, then PDEs have derivatives 
only with respect to ti (but not with respect to x,). We may take F(z , Zi, . . . , Zn) = F(z ) in 
the eq.(228), so that eq.(228) gives us the algebraic equation (41). 

Solutions to Pohlmeyer equation. One must substitute p( mj ) = A5 m j, = 0, N = 2 
into the eq.(228), which yelds the algebraic equation (61) describing solution manifold to the 
eq.(17). Solution manifold to the eq.(20a) corresponds to the reduction 

F(z ,qi,q 2 ) = F(z ,qi)6(q2- z qi), x 2 = 0, (230) 

which yields the algebraic equation (64). 



3.3 Second version of the dressing method for PDEl(^i,^2; w^) with 
reduction (35,37). 

Disadvantage of the dressing algorithm described above is a poor available solutions space to 
PDEOs, since all freedom is assotiated with the function R which is Q x Q matrix function, while 
w has rioQ 2 scalar fields. However, we know that the algebraic approach provides a fullness of 
the solutions space. Thus we may ask if there is another version of the dressing method which 
yelds fullness of the solutions space as well. 
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Such version of the dressing method is based on the integral equation (114) as n Q x noQ 
matrix equation (it was Q x Q matrix equation in the above algorithm). For simplicity, we 
consider only those PDEOs which are assotiated with PDEls. Thus, PDEO (147) with nonzero 

p (m0) ig 

not the subject of this subsection. 
This version of the dressing algorithm is very similar to one discussed in Sees. 3. 1,3. 2. In 
particular, all equations for dressing and spectral functions remain correct, only matrix dimen- 
sionalities must be changed. So we just underline a few features. 

1. It follows from the first paragraph of this section, that PDE0(t m ; n , w ), m — 1,2, with 
any fixed n is assotiated with its own integral equation (114). 

2. Field 

w = H 2 *U*H 1 (231) 

must have the structure (21a). As before, is a solution of PDE1. However, are 
not necessary scalars unlike Sees. 3.1 and 3.2. Since structure (21a) must be supported by 
the nonlinear PDEs (8), parameters should be defined as in the algebraic approach, 
see a paragraph after the eqs.(21). 

3. In order to establish relation between PDEl(£i, t 2 ; w^) and PDE0(t m ; n , w), m = 1,2, 
we take the reduction (141) in the next form: 

A*Hi = H 2 * U * A * Hi = 0. (232) 

Now the functions A, A and A are defined as follows (compare with the formulae (157)): 

A(\,fi) = A(\,fi) = \5(\-fi)I no , A(\,fx) = \5(\-fx)I 2no . (233) 

4. Similarly to the n.3, in order to establish relation between PDEl(ti, t 2 ; w^) and PDE0(t m ; /co- 
rn = 1, 2, we take the reduction (142) in the next form: 

H 2 *A = H 2 * A * U * Hi = 0. (234) 

5. Owing to n.3, derivation of PDE0(t m ; no, w) in this section is equivalent to derivation 
made in Sec. 3. 1.2 with p( m °) = 0, see resulting equation (147). This fact suggests us the 
general form of PDE0(t m ; no, w), eq.(25). 

6. Similarly, owing to n.4, we may derive general form of PDE0(t m ; ko,w), see eq.(26). We 
do not represent appropriate calculations. 

7. Eq.(228) must be replaced by the next equation 

w aP = J2[ e m=1J=1 F ai (w,XiI no -J2p { ^^m,---, (235) 

7=1 m=l 

— v E E T^p^(w)t„ 



,m — 1 j — 1 



m=l 



7/3 



a = l,...,Q, (3 = 1,.. .,n Q. 



We see that the algorithm discussed in this section provides the same richness of the solutions 
space as the algebraic approach does. Namely, we have n^Q 2 arbitrary functions F a p, a = 
1, . . .,Q, (3 = 1, . . .n Q, and diagonal (T^ ^ 0) or arbitrary (T^ = 0)) matrix parameters 
r {j) , j > 1. 
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4 Conclusions 



First of all, we have represented the algorithm for the derivation of the system of non-differential 
equations which implicitly describes the special type of solutions to Self-dual type S-integrable 
equations (11). This is possible due to the relation between S-integrable PDEs (PDEls) and 
appropriate family of matrix first order quasilinear PDEs (PDEOs). 
As we have seen, PDEOs have set of remarkable properties: 

1. Manifold of PDEOs is splitted into subclasses. Each subclass joins such PDEOs which are 
lower dimensional reductions of appropriate PDE1. These PDEOs admit infinitely many 
commuting flows, similar to PDEls. Besides, there is a subclass of PDEOs which is not 
related with PDEls; we have not identified the commuting flows for this subclass. 

2. Any PDEO admits the spectral system consisting of two linear equations for the spectral 
function V(A; x). But only one of these equations has derivatives with respect to variables 
Xi and ti, unlike the spectral systems for PDEls. This fact gives rise to all following 
properties written in this list. 

3. The solutions space of PDEOs is implicitly described by the system of non-differential 
equations. This is a typical description for the first order quasilinear equations [19, 20, 
21, 29, 30]. 

4. As we know, PDEls (like any S'-integrable PDE) may be written using one of the following 
procedures: (a) applying H2 * A n * and *A> * H\ to the spectral problem and using the 
definitions of the dependent variables in terms of the spectral and dressing functions, 
see eq.(140), or (b) as the compatibility conditions of the appropriate spectral problems. 
However, PDEOs may be derived only using the first procedure. 

5. Functions p^\z) and p( m: >\z) may be replaced by p^ m ^(z; x) and p^ m ^(z; x) respectively, 
which are representable by positive power series of z with scalar coefficients explicitly 
depending on x. The solutions spaces to such PDEs is implicitly described by the system 
of the non-differential equations supplemented by the pair of decoupled linear PDEs with 
non-constant coefficients which, in turn, are integrable by the method of characteristics. 
We leave this statement without proof. 

Although all basic conclusions have been established using the algebraic approach, we rep- 
resent a version of the dressing method allowing to derive the same system of non-differential 
equations describing solutions space to PDEOs and PDEls. In addition, this method reproduces 
the manifold of classical solutions to PDEls. 

The fact that the dressing method exhibits a connection with the method of characteristics 
is important for unification of the methods for solving of nonlinear PDEs in multidimensions. 
It tells us that the dressing method is more flexible in comparison, at least, with some other 
integration methods mentioned in this paper, such as Inverse Spectral Transform (1ST) (S- 
integrability) , direct linearization (C-integrability) and method of characteristics. In fact, direct 
linearization may not be applied to both S'-integrable equations and PDEOs. Similarly, 1ST may 
not be applied to PDEOs and to C-integrable equations. Method of characteristics was modified 
to hodograph and generalized hodograph method allowing to study some (S-integrable) PDEs 
by the method of hydrodynamic reductions. However, using this technique we face a problem 
of the explicit representation of solutions. 
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On the contrary, although originally the dressing method was invented as a special tech- 
nique for solving S'-integrable equations, it was shown in [30] as well as in this paper that the 
properly modified dressing method works also for PDEOs. It yelds exactly the same matrix 
non-differential equation describing solutions space as method of characteristics does. The fact 
that the dressing method may be a starting point for solving C-integrable equations has been 
demonstrated in [54] as well as in [30]. 

Author thanks Professors P.M.Santini and A.B.Shabat for useful discussions and Professor 
I.M.Krichever for important comments. The work was supported by INTAS Young Scientists 
Fellowship Nr. 04-83-2983, RFBR grants 04-01-00508, 06-01-90840, 06-01-92053 and grant NS 
7550.2006.2. 
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